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Abstract  of  Dissertation  Presented  to  the  Graduate  School 
of  the  University  of  Florida  in  Partial  Fulfillment  of  the 
Requirements  for  the  Degree  of  Doctor  of  Philosophy 

ANALYTIC  ONE-ELECTRON  RESPONSE  PROPERTIES, 

MOLECULAR  GRADIENTS,  AND  FORCE  CONSTANTS 
IN  MANY-BODY  METHODS 

By 

Edward  Alan  Salter 
August,  1988 

Chairman:  Rodney  J.  Bartlett 

Major  Department:  Chemistry 

The  theory  of  analytic  first  and  second  derivatives  of  the 
Coupled-Cluster  (CC)  energy  is  presented.  First  .derivatives  of  the 
energy,  corresponding  to  first-order  one-electron  response  properties 
and  molecular  gradients,  are  expressed  in  terms  of  the  response  (or 
"relaxed")  density  for  efficient  evaluation.  Explicit  expressions  for 
the  derivative  energy  and  response  density  for  the  Coupled-Cluster 
Singles,  Doubles  and  Triples  (CCSDT)  model  are  presented;  expressions 
for  the  finite-order  Many-Body  Perturbation  Theory  (MBPT)  models 
through  full  fourth  order  are  extracted  from  low-order  iterations  of 
the  derivative  CCSDT  expressions. 

Second  derivatives  of  the  CC  energy,  corresponding  to  second-order 
one-electron  response  properties  and  force  constants,  are  expressed  in 
terms  of  the  response  density  and  first-derivative  CC  t-amplitudes . 
Explicit  expressions  for  the  second  derivative  of  the  energy  are 
presented  for  the  CCSDT-1  model,  the  CCSDT  model  truncated  to  linear 
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triple-excitation  contributions.  Expressions  for  the  second 
derivative  of  the  finite-order  MBPT  energies  through  full  fourth  order 
are  presented  as  well. 
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BACKGROUND:  MODELS  FOR  CORRELATION 


Many  of  the  traditional  tools  for  approximating  eigenfunctions  and 
eigenvalues  of  the  electronic  Hamiltonian  are  based  upon  the 
variational  principle.  Given  the  choice  of  an  atomic  orbital  basis 
set,  which  usually  consists  of  Gaussian  functions  centered  upon  the 
nuclei  in  the  molecular  system,  the  Self-Consistent  Field  (SCF)  method 
is  commonly  applied.  The  SCF  procedure  determines  a set  of  molecular 
orbitals,  one-electron  functions  of  space  and  spin,  as  linear 
combinations  of  the  atomic  orbital  basis  functions.  The  SCF 
wavefunction  is  an  antisymmetrized  product  (Slater  determinant)  of  the 
occupied  molecular  orbitals;  the  SCF  electronic  energy  is  the 
expectation  value  of  the  electronic  Hamiltonian  over  the  SCF 
wavefunction.  The  SCF  molecular  orbital  coefficients  are 
variationally  optimal  for  a single-determinant  description  of  the 
electronic  wavefunction.  In  the  SCF  model,  each  electron  moves  in  an 
effective  potential — an  average  of  its  interactions  with  all  other 
electrons  in  the  system,  plus  its  coulomb  attraction  for  the  nuclei. 
The  SCF  method  is  inexpensive  but  not  very  accurate  in  its  prediction 
of  molecular  structure  and  properties.  Consequently,  there  exist 
"correlated"  methods  which  restore  the  true  electronic  interactions  by 
going  beyond  the  single-determinant  description  of  the  electronic 
wavefunction  and  the  effective  SCF  Hamiltonian. 

As  a by-product  of  the  SCF  method,  molecular  orbitals  of  two  types 
are  determined:  the  "occupied"  orbitals  and  the  "virtual"  (or 
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"unoccupied")  orbitals.  The  virtual  SCF  orbitals,  in  conjunction  with 
the  SCF  determinant,  are  used  to  define  "excited"  determinants.  The 
collection  of  all  determinants  obtained  by  replacing  one  occupied 
orbital  in  the  SCF  determinant  with  a virtual  orbital  compose  the 
singly-excited  determinants.  The  collection  of  determinants  where  two 
occupied  orbitals  have  been  replaced  by  two  virtual  orbitals  compose 
the  doubly-excited  determinants.  Triple  and  higher  excitations  are 
analogously  defined.  The  Configuration  Interaction  (Cl)  method  refers 
to  the  construction  and  subsequent  diagonalization  of  the  matrix 
representation  of  the  Hamiltonian  in  the  basis  of  all  such 
determinants.  The  Cl  electronic  wavef unction  is  then  the  SCF 
determinant  plus  a linear  combination  of  excited  determinants. 
Furthermore,  the  Cl  wavefunction  is  variationally  the  best  possible 
wavef unction  in  the  basis  set  spanned  by  the  determinants,  and  by 
definition,  is  "fully  correlated." 

For  practical  purposes,  the  full  Cl  problem  is  not  tractable 
because  the  complete  set  of  excited  determinants  is  usually  extremely 
large.  Indeed,  if  the  number  of  atomic  basis  functions  (M)  is  large, 
there  are  approximately  M!  excited  determinants.  The  basis  set  of  the 
Hamiltonian  is  therefore  truncated  to  include  only  excited 
determinants  of  certain  types.  For  example,  simple  models  might 
contain  doubly-excited  determinants  (CID)  or  singly-  and  doubly- 
excited  determinants  (CISD).  The  lowest  eigenvalue  and  corresponding 
eigenvector  are  extracted  from  the  truncated  Cl  Hamiltonian  matrix, 
yielding  a wavefunction  that  is  variationally  the  best  wavefunction  in 
the  truncated  space. 

Alternatively,  the  deficiencies  of  the  SCF  wavefunction  and  SCF 
effective  Hamiltonian  are  circumvented  by  the  Multi-Configuration  SCF 
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(MCSCF)  method.  Like  the  truncated  Cl  vavefunction,  the  MCSCF 
vavefunction  is  a linear-combination  of  a selected  set  of 
determinants,  where  the  coefficients  are  variationally  optimal. 

Unlike  the  Cl  method,  however,  the  orbitals  are  not  fixed  by  the  SCF 
method  but  are  allowed  to  vary  as  well. 

Many-body  methods  such  as  the  Coupled-Cluster  (CC)  method  and 
Many-Body  Perturbation  Theory  (MBPT)  are  often  the  modern  methods  of 
choice  for  solving  the  correlation  problem.  Although  many-body 
methods  are  not  variational,  they  provide  accurate  approximations  to 
the  full  Cl  wavefunction  and  full  Cl  energy  eigenvalue  and  possess  the 
size-extensive  property  of  the  full  Cl  model.  A model  is  said  to  be 
size-extensive  if  the  energy  scales  properly  with  the  size  of  the 
molecule.  Consider,  for  example,  a system  composed  of  N identical, 
non-interacting  components.  The  full  Cl  energy  of  the  system  will 
necessarily  be  N times  the  energy  of  a single  component.  The 
approximate  many-body  energy  of  the  system  will  also  be  N times  the 
approximate  many-body  energy  of  a single  component.  Chemical  systems 
are  understood  at  the  most  fundamental  level  of  approximation  to  be 
composed  largely  of  independent  electron  pairs;  thus,  the  errors  of 
truncated  Cl  or  MCSCF  approximations,  which  do  not  preserve  the  size- 
extensive  property,  can  be  quite  large. 


INTRODUCTION 


The  evaluation  of  energy  derivatives  is  a necessary  element  of  the 
ab  initio  prediction  of  molecular  structure  and  properties. 

Derivatives  of  the  energy  with  respect  to  nuclear  displacement 
parameters,  called  "molecular  gradients",  represent  the  forces  acting 
on  the  nuclei  and  are  necessary  input  for  search  algorithms  which 
locate  critical  points  on  a molecule's  potential  energy  surface.  Once 
a critical  point  has  been  found,  the  second  derivatives  of  the  energy, 
evaluated  perhaps  by  finite-difference  of  analytic  molecular 
gradients,  characterize  the  point  as  a minimum  (a  stable  molecular 
structure)  or  a saddle  point  (a  transition  state)  or  neither.  Second 
derivatives  of  the  energy  make  possible  the  prediction  of  harmonic 
vibrational  frequencies  and  intensities,  as  well.  Derivatives  of  the 
energy  with  respect  to  external  field  perturbation  parameters  are 
directly  associated  with  the  one-electron  response  properties  of  the 
molecule.  The  prediction  of  molecular  structure  and  properties  using 
highly  correlated  many-body  methods  have  become  standard  in  quantum 
chemical  studies;  in  this  work,  procedures  for  the  efficient 
evaluation  of  analytic  first  derivatives  of  the  MBPT/CC  energy  are 
presented. 

One-electron  response  properties  are  often  computed  in  the 
framework  of  approximate  Cl  models  as  derivatives  rather  than  as 
expectation  values.  In  the  full  Cl  model,  the  derivative  approach 
yields  the  same  result  as  an  expectation  value;  a Hellman-Feynman 
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Theorem  is  said  to  exist  for  the  full  Cl  wavefunction.  The  Hellman- 
Feynman  Theorem  applies  to  the  SCF  model,  too.  For  many  approximate 
models  of  correlation,  however,  the  expectation  value  and  derivative 
are  not  the  same.  In  fact,  the  derivative  is  the  expectation  value 
plus  additional  terms  which  usually  bring  the  derivative  closer  to  the 
full  Cl  property  value.  First-order  response  properties  in  the 
MBPT/CC  model  are  shown  to  be  simple  dot-products  between  the  response 
or  "relaxed"  MBPT/CC  density  matrix  and  the  one-electron  perturbation 
matrix  of  interest. 

Molecular  gradients  are  treated  as  special  cases  of  a response 
property,  wherein  the  basis  functions  are  dependent  upon  the  nuclear 
displacement  perturbation.  The  Gaussian  basis  functions  are  fixed  to 
the  nuclei,  and  when  the  nuclei  move,  the  basis  functions  move  as 
well.  In  addition  to  the  response-property  term,  there  are  terms 
involving  the  derivative  of  the  atomic  orbital  basis  functions  in  two- 
electron  integrals  and  in  overlap  integrals. 

The  evaluation  of  energy  derivatives  of  variational  methods  can  be 
simplified  in  ways,  which,  in  general,  the  evaluation  of  energy 
derivatives  of  many-body  methods  cannot.  The  energy  is  generally  a 
function  of  the  molecular  orbital  coefficients  c,  the  coefficients  in 
the  determinantal  expansion  of  the  wavefunction  C,  and  the  one-  and 
two-electron  integrals  (I)  involving  the  atomic  orbital  basis  set 
functions.  The  derivative  of  the  energy  is  in  general: 

3E(c(X) ,C(X) » I(X) ) 9E  3c  3E  3C  3E  31 

' = + ' + • 

3X  9c  3X  9C  3X  91  9X 

In  the  case  of  the  SCF  model,  the  wavefunction  is  a single  determinant 
and  the  energy  is  optimal  with  respect  to  changes  in  the  molecular 
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orbital  coefficients;  that  is,  3E/3c  = 0.  It  is  not  necessary  to 
evaluate  3c/3X,  the  derivative  of  the  molecular  orbital  coefficients 
with  respect  to  the  perturbation  parameter.  If  X represents  a nuclear 
displacement  parameter,  the  third  term  indicates  that  it  is  necessary 
to  evaluate  3E/3I,  how  the  SCF  energy  varies  relative  to  the  one-  and 
two-electron  integrals,  and  in  turn  to  evaluate  31/ 3X,  how  the 
integrals  vary  relative  to  the  nuclear  displacement  parameter.  Such 
contributions  distinguish  the  Generalized  Hellman-Feynman  Theorem  from 
from  the  ordinary  Hellman-Feynman  Theorem.  Similarly,  in  the  case  of 
the  MCSCF  model,  we  have  3E/3C  = 0 and  3E/3c  =0,  and  again  only  the 
third  term  remains.  In  Cl  methods,  the  energy  is  optimal  with  respect 
to  the  determinantal  coefficients  only;  that  is,  we  have  3E/3C  - 0 and 
3E/3c  * 0. 

In  the  Coupled-Cluster  method1-11  and  its  finite-order  MBPT 
approximations , 5 , 8 , 12  14  energy  is  not  optimal  with  respect  to 
the  molecular  orbital  coefficients  or  the  determinantal  coefficients. 
It  is  necessary  to  take  into  account  the  derivative  of  the  CC  energy 
with  respect  to  the  SCF  coefficients  and  with  respect  to  the 
determinantal  coefficients,  and  in  turn,  the  derivative  of  the  SCF 
coefficients  and  the  determinantal  coefficients  with  respect  to  the 
perturbation  parameter.  It  might  appear  that  for  the  evaluation  of 
several  properties  and/or  the  3N  molecular  gradients  (for  a system  of 
N nuclei),  it  is  necessary  to  evaluate  3C/3X  and  3c/3X  for  each 
perturbation  of  interest.  However,  the  equation  governing  the 
response  of  the  SCF  coefficients  is  a linear  equation,  as  is  the 
equation  governing  the  response  of  the  CC  determinantal  coefficients. 
Consequently,  the  first  and  second  terms  of  the  energy  derivative 
expression  can  be  transformed  into  a computationally  efficient  form. 
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The  determinantal  coefficient  derivatives  for  each  perturbation  need 

not  be  determined;  only  a single,  perturbation-independent  linear 

equation  must  be  solved.  Analogously,  the  SCF  coefficient  derivatives 

for  each  perturbation  need  not  be  determined;  once  again,  it  is 

necessary  to  solve  only  a single,  perturbation-independent  linear 

equation  in  order  to  evaluate  the  energy  derivative.  A special 

property,  the  response  or  "relaxed"  density,  emerges  naturally  from 

the  transformation  to  perturbation-independent  expressions. 

Analytic  gradients  for  the  SCF,1"^  MCSCF,1^  Cl,12  and  MBPT(2) 

models  have  been  available  for  some  time.  Recent  developments  in  the 

analytic  evaluation  of  CC/MBPT  gradients  have  been  based  upon 

perturbed  CC  theory.19-22  The  first  reports  of  the  theory  and 

23  24 

application  of  MBPT(3),  D-MBPT(4)  and  CCD  gradients  ’ were  not 

suitable  for  routine  applications  since  they  did  not  take  advantage  of 

the  linear  nature  of  the  CC  derivative  equations.  The  general  theory 

25  26 

of  CC  derivatives  was  later  presented  in  an  efficient  form,  ’ and 

27 

analytic  CCSD  gradients  were  subsequently  reported.  Analytic  SDQ- 

28 

MBPT(4)  gradients  have  been  recently  reported,  but  the 

implementation  did  not  take  advantage  of  the  linear  nature  of  the  SCF 

27  28 

derivative  equations.  The  most  recent  gradient  work  ’ has 

29 

incorporated  the  discussion  of  Handy  et  al.  on  the  elimination  of 
singularities  which  arise  when  there  are  degenerate  or  near-degenerate 
orbital  eigenvalues. 

Since  the  equations  governing  the  second-order  response  of  the  SCF 
coefficients  and  the  CC  determinantal  coefficients  are  linear,  there 
is  no  need  to  solve  explicitly  for  second  derivative  SCF  coefficients 
or  second  derivative  determinantal  coefficients  in  order  to  evaluate  a 
CC  second  derivative.  In  simplest  form,  second  derivatives  of  the  CC 
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energy  are  expressed  in  terms  of  the  relaxed  density  and  first- 
derivative  CC  t-amplitudes . 

18 

The  evaluation  of  second  derivatives  at  the  SCF  level  of  theory 

is  commonplace  today;  the  theory  and  application  of  analytic  many-body 

second  derivatives  is  currently  limited  to  the  simplest  model — second 

29  30 

order  perturbation  theory.  ’ 

In  this  work,  the  CCSDT  derivative  energy  expressions  are  derived 
using  diagrammatic  rules  developed  for  second-quantized  operators. 

The  finite-order  MBPT  derivative  models  through  full  fourth  order  are 
obtained  by  iteration  of  the  CC  derivative  equations.  The  theory  of 
CC/MBPT  second  derivatives  is  presented.  Explicit  expressions  for  the 
second  derivative  of  the  energy  are  presented  for  the  CCSDT-1  model, 
the  CCSDT  model  truncated  to  linear  triple-excitation  contributions. 
Second  derivative  expressions  for  the  finite-order  MBPT  models  through 
full  fourth  order  are  provided  as  well. 


THE  COUPLED-CLUSTER  MODEL 


Consider  the  Schrbdinger  Equation  corresponding  to  a molecular 
system,  where  the  electronic  Hamiltonian  H(X)»  its  eigenfunction  T(X)> 
and  its  eigenvalue  E(X)  are  dependent  upon  some  parameter  X : 

H(X)  Y(X)  = E(X)  Y(X)  (1) 

2 -1  -1 
H(X>  - -1/2  I 7,  - I Z„rol  . £ XQ,  * £ r • 

i i,a  i 

The  parameter  may  control  the  displacement  of  an  atomic  nucleus  in  the 
molecular  system,  or  it  may  be  the  strength  of  an  external  one- 
electron  potential  energy  perturbation  Q^,  such  as  a uniform  electric 
field  applied  in  the  z-direction:  XQi  = Xzr  The  one-body  part  of  the 
electronic  Hamiltonian  is  composed  of  three  parts:  the  kinetic  energy 
operators  of  the  electrons;  the  potential  energy  operators  of  the 
electrons  due  to  the  nuclei  present;  and  a potential  energy  operator 
due  to  an  externally  applied  field,  if  present.  The  two-body  operator 
is  the  electron-electron  repulsion  operator.  Here,  the  labels  i and  j 
represent  the  electrons  in  the  system;  the  label  a represents  the 
nuclei,  and  is  the  charge  of  nucleus  a.  In  correlated  methods,  we 
seek  the  solution  to  Equation  (1)  in  the  space  of  the  N-electron 
determinants  (antisymmetrized  products)  formed  from  a basis  of 
orthonormal  one-electron  spin-orbitals  {<KX)}»  The  solution  is  built 
upon  an  approximate  single-determinant  wavefunction,  called  the 
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reference  function  and  denoted  by  |0>.  The  total  energy  of  the 
molecular  system,  within  the  Born-Oppenheimer  or  "fixed  nuclei" 
approximation,  is  obtained  by  adding  the  nuclear  repulsion  energy  to 
the  electronic  energy  E(X)- 

The  reference  energy  is  defined  as  the  expectation  value  of  the 
Hamiltonian  over  the  reference  vavefunction  <0|H(X)|0>.  Subtraction 
of  the  reference  energy  from  each  side  of  Equation  (1)  yields: 

(H(X)  - <0 | H ( X) | 0> ) Y(X)  = (E(X)  - <0|H(X)|0>)  T(X) 

or  Hn(X)  Y(X)  = AE(X)  *(X)  • (2) 

H^(X)  is  the  "normal  product"  form  of  the  Hamiltonian  and  is 
represented  in  terms  of  the  second-quantized  operators  by: 


Hn(X)  = l 

p.q 


h + Z <pi | | qi> 

pq 


{p+q}  + - Z <pq||rs>{p  q sr}  (3) 
4 

p»q. 


r , s 


By  convention,  the  indices  i,j,k,l,...  represent  the  spin  orbitals 
which  are  occupied  in  the  reference  function  |0>.  The  indices 
a,b,c,d,...  correspond  to  the  virtual  (or  unoccupied)  spin  orbitals. 
The  indices  p,q,r,s,...  are  unrestricted. 

The  symbols  hpq  and  <pq||rs>  in  Equation  (3)  represent  integrals 
wherein  the  one-  and  two-electron  operators  of  the  Hamiltonian  act 
upon  the  orbital  basis  functions: 


' ★ 

' 2 -1 

h 

pq 

v1) 

-1/2  V1  - Z Z r . + XQ-, 
1 a al  1 

l a 

( 1 ) 3ti  3 o-, 
q v ' 1 1 
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<pq||rs>  = 


y1*  Y2) 


r12  P12 


4>r ( 1 ) 4>s(2)  0T13ff13T23ff2 


The  variables  and  Tj  are  the  coordinates  of  spin  and  space  of 
electron  number  one,  respectively.  The  operator  P^2  is  the 
permutation  operator;  thus  the  integral  <pq||rs>  is  antisymmetric  with 
respect  to  the  exchange  of  p with  q or  r with  s.  The  integral  h is 
explicitly  dependent  upon  X if  a one-electron  potential  energy 
perturbation  is  present.  Of  course,  if  the  perturbation  is  a 
nuclear  displacement,  the  h integrals  depend  upon  X through  the 
operator  Z^r^  as  well.  In  general,  both  h and  <pq||rs>  are 
implicitly  dependent  upon  X through  the  orbital  basis  functions,  since 
the  basis  functions  vary  with  molecular  geometry  or  external  field 
strength.  Consider  the  common  choice  of  SCF  molecular  orbitals  as  the 
orbital  basis.  Self-Consistent  Field  orbitals  are  determined  in  the 
presence  of  the  effective  Fock  potential  which  contains  the  entire 
one-electron  part  of  the  Hamiltonian.  In  addition,  the  SCF  molecular 
orbitals  are  determined  as  linear  combinations  of  primitive  AO  basis 
functions,  and  if  the  perturbation  is  a nuclear  displacement,  the 
primitive  AO  basis  functions  centered  upon  that  nucleus  are  displaced 
as  well. 

The  operators  p+  and  q are  the  second-quantized  creation  and 
annihilation  operators,  respectively.  Excited  determinants  are 
represented  in  second-quantization  by  a string  of  creation  and 
annihilation  operators  acting  on  the  reference  function,  whereby  an 
electron  in  an  occupied  orbital  is  "annihilated"  and  an  electron  is 
subsequently  "created"  in  a virtual  orbital: 
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{a+i}|0>  = 

|S> 

Singly-excited 

determinant 

{a+ib+j } | 0>  = 

1 D> 

Doubly-excited 

determinant 

{a+ib+jc+k} |0>  = 

| T> 

Triply-excited 

determinant 

etc. 

Brackets  enclosing  a second-quantized  operator  string  {ABCD...} 
denote  a normal  product  of  operators.  The  normal  product  of  a second- 
quantized  operator  string  is  defined  by  moving  all  a,b,...  and 
i+,j+,...  operators  to  the  right  of  the  operator  string  by  use  of  the 
anticommutation  relations  which  govern  second-quantized  operators: 

[a+,b]+  = Sab  and  [i+,j]+  = S.j  . 

The  anticommutation  relations  follow  from  the  defining  conditions: 

31 

a|0>  = i+|0>  = 0.  The  normal-ordered  form  of  the  Hamiltonian  in 

Equation  (3)  is  obtained  by  application  of  Wick's  Theorem  to  the 

ordinary  second-quantized  Hamiltonian.  Wick's  Theorem  states  that  a 

string  of  m second-quantized  operators  ABCD...  equals  the  normal 

product  of  the  string  plus  all  possible  one-fold,  two-fold,...,  up  to 

31 

m-fold  contracted  normal  products  of  the  string. 


ABCD. . . 


n n n 

= {ABCD...}  + {ABCD...}  + {ABCD...}  + {ABCD...} 

+ (all  other  one-fold  contracted 


+ {ABCD...} 
products) 


n n 
+ {ABCD 


,}  + {ABCD. 


n 


47 


} + {ABCD...}  + {ABCD...} 
+ (all  other  two-fold  contracted  products) 


+ (all  m-fold  contracted  products) 

A contraction  between  operators  A and  B is  defined  as 

AB  = AB  - {AB}  . 

The  only  nonvanishing  contractions  are 

m n + 

i j = 6.  . and  ab  = 5 , . 

J 13  ab 

It  follows  that  all  contractions  of  a normal  product  of  operators 
{ABCD...}  are  vanishing  ones.  The  expectation  value  of  a second- 
quantized  operator  string  over  the  reference  (or  Fermi  vacuum)  can  be 
evaluated  by  forming  all  of  the  fully  contracted  products. 

<0| ABCD. . . |0>  = (all  m-fold  contracted  products) 

So,  the  expectation  value  of  a normal  product  over  the  reference  is 
zero. 


<0 | {ABCD. . .} |0> 


0 


14 


31 

The  Generalized  Wick's  Theorem  states  that  a product  of  operator 
strings  in  which  some  operator  strings  are  already  in  normal  product 
form  is  given  as  the  normal  product  of  the  overall  string  plus  the  sum 
of  all  possible  normal  products  with  contractions;  however,  no 
contractions  between  pairs  of  operators  within  an  original  normal 
product  are  included.  A corollary  is  that  for  the  expectation  value 
of  a product  of  normal  product  operators  over  the  reference  (or  Fermi 
vacuum),  only  the  fully  contracted  products  are  non-zero. 

The  representation  of  the  Fock  operator  in  the  orbital  basis  is 
given  by  the  matrix  elements 


f 

pq 


h + E <pi | | qi>  • 

pq 

i 


So  we  may  represent  HN(X)  as 


Hn(X)  = E fpq{p+q)  + - 2 <Pq I I rs> {p  q sr} 


1 

4 


p,q 


p»q» 

r , s 


(4) 


or  Hn(X)  = fN(X)  + VN(X)  • 

If  the  SCF  wavefunction  is  chosen  as  the  reference  function, 

<0 | H( X) I 0>  is  the  SCF  energy,  AE(X)  is  the  correlation  energy,  and  the 
Hamiltonian  simplfies  to 

Hm(X)  = Z 2 f hU+b)  + " L <Pql  |rs>{p+q+sr)  , 

N ij  aD  4 

i,j  a,b  p,q, 

r , s 
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since  the  matrix  elements  f j and  f^a  are  zero  for  SCF  orbitals.  For 

the  particular  choice  of  canonical  SCF  orbitals,  we  have,  by 

definition,  f ^ and  f^  = $abea»  where  ep  denotes  the 

molecular  orbital  eigenvalue  of  orbital  p. 

1 —1  1 

The  CC  method  procedes  by  choosing  an  exponential  form  of 
Y(X).  That  is, 


V(X)  = eT(X) |0>  , 


(5) 


where  T(X)  is  a sum  of  excitation  operators 


T(X)  = TX(X)  + T2(X)  + T3(X)  + 


The  excitation  operators  are  given  by 


T, (X)  = E t (X)  {a+i} 

1 i 

i,a 


t2(x) 


Tn(X) 


1 

4 


ab 

: i 

ij 


- I t (X)  {a+ib+j} 


a,b 


(n!)- 


ab ... 

ij 


Z t (X)  {a+ib+j . . . } 


i.l » 

a,b, 


where  the  t-amplitudes  are  to  be  determined. 


Substitution  of  Equation  (5)  into  Equation  (2)  yields: 
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Hn(X)  eT(X)|0>  = AE(X)  eT(X)|0> 


or 


e"T(X)HN(X)  eT(X)|0>  = AE(X)|0>  . 


(6) 


-A  A 

The  Hausdorff  expansion  of  e Be  yields  a series  of  commutators: 


1 


1 


e~ABeA  = B - [A,B]  + - [A, [A,B] ] - - [A, [A, [A, B] ] ] + ...  . 

2!  3! 

Since  a product  of  normal  product  operators  can  be  written  as  a sum  of 
contracted  terms,  it  follows  that  a product  AB  has  a "connected"  part 
and  a "disconnected"  part.  In  the  connected  part  there  is  at  least 
one  contraction  that  connects  A to  B;  in  the  disconnected  part  there 
is  no  such  contraction.  For  a commutator  of  operators  [A,B],  the 
disconnected  parts  cancel,  leaving  only  the  connected  parts: 


[A,B]  = AB  - BA  = (AB)c  - (BA)c  + (AB)d  - (BA>d  = [A,BJc  . 


Thus  the  expansion  becomes 


Furthermore,  if  A is  an  excitation  operator,  A cannot  give  rise  to  a 
connected  contribution  when  it  is  to  the  left  of  other  operators. 
Consequently,  the  expansion  reduces  to  the  following  form: 
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e"ABeA  = B + (BA)  + - ((BA)  A)  + - (((BA)  A)  A)  + ... 

c 2!  c c 3!  c c c 

= (Be\  • 

So  Equation  (6)  can  be  represented  as 

(Hn(X)  eT(X))c|0>  = AE(X)|0>  . (7) 

By  projecting  Equation  (7)  on  the  left  by  <0 | , we  obtain  the  CC 
energy  expression: 


AE(X)  = <0|(Hn(X)  eT(X))c|0>  . (8) 

The  t-amplitudes  are  determined  by  solving  the  set  of  coupled 
equations  obtained  by  projection  of  the  singly-,  doubly-, triply-  ... 
excited  determinants  (S,  D,  T,  ...): 

<S|(Hn(X)  eT(X))c|0>  = 0 (9) 

<D|(Hn(X)  eT(X))c|0>  = 0 

<T|(Hn(X)  eT(X))c|0>  = 0 


Various  CC  models  result  from  choices  of  truncation  of  T(X)>  For 
example,  the  CCD  model  is  defined  by  the  choice  T(X)  = T2(X)»  the  CCSD 
model  is  defined  by  T(X)  = T^(X)  + T2(X)  and  the  CCSDT  model  is 
defined  by  T(X)  = T^X)  + T2(X)  + T3(X).  In  the  limit  of  including 
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all  excitation  operators  in  T(X),  the  full  Cl  vavefunction  and  energy 
is  recovered. 

Diagrammatic  rules  are  used  to  identify  all  of  the  fully 
contracted  products  in  the  CC  amplitude  equations  and  energy 
expression;  the  diagrams  are  then  interpreted  algebraically.  The 
operators  fN(X)  and  WN(X),  which  compose  the  normal  product 
Hamiltonian  HN(X),  are  represented  diagrammatically  in  Figure  (1). 

The  CC  excitation  operators  T^X),  T2(X),  and  T3(X)  are  represented  by 
the  diagrams  in  Figure  (2).  The  diagrams  representing  each  of  the 
fully  contracted  products  in  the  CC  energy  expression  is  shown  in 
Figure  (3).  The  diagrammatic  equations  for  the  T^(X)>  T^(X)  and  T^(X) 
amplitudes  are  shown  in  Figures  (4),  (5),  and  (6),  respectively.  In 
the  diagrams,  the  operator  diagrams  from  H^(X)  and  e are  connected 
together;  each  connection  represents  a contraction  between  the 
operators.  A diagram's  excitation  level— the  number  of  pairs  of 
unconnected  lines  at  the  top  of  the  diagram,  matches  the  excitation 
level  of  the  projected  determinant.  That  is,  all  diagrams  in  the 
energy  expression,  obtained  by  left-projection  of  <0|,  have  a zero 
excitation  level.  The  diagrams  in  the  T^X)  equation,  obtained  by 
left-projection  of  | S>  have  an  excitation  level  of  one.  The  T2(X) 
equation  consists  of  diagrams  of  excitation  level  two,  and  so  forth. 
The  explicit  algebraic  expressions  for  the  CCSDT  model  are  presented 
elsewhere. ^ 

For  the  choice  of  canonical  SCF  orbitals  as  the  orbital  basis,  the 
terms  involving  the  fN(X)  operator  simplify  because  f is  zero  except 
along  the  diagonal  where  the  orbital  eigenvalues  are.  Consider  an 
iterative  form  of  the  equations  where  the  diagrams  (3a, b)  in  Figure 
(4)  are  taken  to  the  other  side  of  the  T^(X)  equation,  diagrams  (2a, b) 
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in  Figure  (5)  are  taken  to  the  other  side  of  the  Tj(X)  equation,  and 
diagrams  (3a, b)  in  Figure  (6)  are  taken  to  the  other  side  of  the  T^(X) 
equation.  Solutions  to  the  equations  may  be  constructed  by  forming 
successive  iterates  of  increasing  order  in  the  correlation 
perturbation  WN(X)-  Iteration  of  the  amplitude  equations  and 
subsequent  substitution  into  the  energy  expression  yield  the  finite- 
order  MBPT  energies.  Attention  is  centered  upon  the  T2(X)  amplitudes, 
since  through  fourth-order  MBPT,  only  the  T2(X)  amplitudes  contribute 
directly  to  the  correlation  energy  by  way  of  the  third  diagram  in 
Figure  (3). 

The  first-order  iterate  of  T2(X)  arises  from  the  constant  term  in 
the  T2<X)  equation  represented  by  diagram  (1)  in  Figure  (5).  The 
amplitudes  are  obtained  by  dividing  the  constant  term  by  the 
eigenvalues  derived  from  the  diagrams  (2a, b),  mentioned  above. 
Substitution  into  the  energy  expression  yields  the  second-order  energy 
MBPT(2) . MBPT  diagrams  are  described  in  terms  of  the  order  of  the 
correlation  perturbation  W^(x)  and  in  terms  of  the  excitation  level  of 
the  determinant  from  which  the  energy  contribution  arises.  The  SCF 
energy  is  complete  through  first  order  in  WN(X);  thus  the  MBPT 
correlation  models  begin  at  second  order  in  U^(X).  The  MBPT(2)  energy 
is  represented  by  a single  MBPT  diagram,  and  it  originates  from 
doubly-excited  determinants  in  the  correlated  wavefunction. 

The  second-order  iterate  of  T2(X)  arises  from  the  substitution  of 
the  first-order  iterate  into  the  linear  terms  of  the  T2(X)  equation- 
-diagrams  (2c-e)  in  Figure  (5).  The  second-order  iterate  is  completed 
by  dividing  through  by  the  eigenvalues  derived  from  the  diagram 
(2a, b).  Substitution  of  the  second-order  iterate  of  T2(X)  into  the 
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energy  expression  yields  the  three  diagrams  of  the  MBPT(3)  model,  all 
arising  from  doubly-excited  determinants,  as  well. 

The  third-order  iterate  of  T2(X)  has  contributions  from  four 

sources:  the  lowest-order  iterate  of  T^(X)>  which  happens  to  be  second 

order  in  VN(X),  contributes  to  T2(X)  by  way  of  diagrams  (3a, b)  in 

Figure  (5);  the  second-order  iterate  of  T2(X)  contributes  to  T2(X) 

through  the  linear  diagrams  (2c-e)  of  Figure  (5);  the  lowest-order 

iterate  of  T.j(X)»  which  also  happens  to  be  second-order  in  VN(X), 

contributes  to  T2(X)  by  way  of  diagrams  (4b, c)  in  Figure  (5);  the 

first-order  iterate  contributes  to  T2(X)  by  way  of  the  quadratic  part 

of  the  T2(X)  equation — diagrams  (5a-d)  in  Figure  (5).  The  result  of 

substitution  of  the  third-order  T2(X)  iterate  into  the  energy 

expression  is  that  corresponding  fourth-order  MBPT  energy 

contributions:  S(4),  D(4),  T(4)  and  Q(4).^  The  letters  S,  D,  T and  Q 

indicate  that  the  contributions  are  generated  by  single,  double, 

triple  and  quadruple  excitations,  respectively.  The  lowest-order 

iterate  of  T^(X)  is  obtained  by  substitution  of  the  first-order 

iterate  of  T2(X)  into  diagrams  (2b, c)  in  Figure  (4)  and  subsequently 

dividing  by  the  eigenvalues  derived  from  diagrams  (3a, b).  The  lowest- 

order  iterate  of  T^(X)  is  obtained  by  substitution  of  the  first-order 

iterate  of  T2(X)  into  diagrams  (la,b)  in  Figure  (6),  and  subsequently 

dividing  by  the  eigenvalues  derived  from  diagrams  (3a, b).  Finally,  we 

note  that  the  quadratic  part  of  the  T2(X)  equation,  giving  rise  to  the 

2 

Q(4)  MBPT  energy  diagrams,  has  its  origins  in  the  T2(X)  part  of  the 
exponential  wavefunction;  the  lowest-order  energy  contributions  from 
quadruply-excited  determinants  are  actually  products  of  doubly-excited 


determinants . 
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For  practical  molecular  systems,  computational  limitations  often 
force  the  restriction  to  a partial  inclusion  of  T^X)  in  the  CC 
equations.  The  CCSDT-1  model  is  defined  by  truncation  of  the  T^(X) 
amplitude  equation  to  diagrams  (la,b)  in  Figure  (6)  and  by  the 
restriction  of  T3(X)  contributions  to  T-^X)  and  T2(X)  to  the  linear 
terms.  That  is,  T^(X)  is  permitted  to  contribute  to  T^(X)  only 
through  diagram  (4)  in  Figure  (4)  and  to  contribute  to  T2(X)  only 
through  diagrams  (4a-c)  in  Figure  (5).  The  lowest-order  T^(X) 
approximation  is  preserved  so  that  the  CCSDT-1  model  is  complete 
through  the  fourth-order  MBPT  energy  and  the  second-order  MBPT 
vavefunction. 

Because  the  equations  for  the  CC  amplitudes  are  connected,  the 
subsequent  substitution  of  the  amplitudes  into  the  energy  expression, 
which  is  connected  and  closed,  yields  only  "linked"  energy  diagrams. 
That  is,  no  diagram  has  a closed,  disconnected  part.  The  size- 

O 

extensive  property  of  the  CC/MBPT  models  follows  from  the  fact  that 
all  energy  contributions  are  linked. 

The  choice  of  a particular  CC  or  MBPT  model  defines  an  approximate 
solution  to  Equation  (2).  The  approximate  total  energy  of  the 
molecular  system  is  obtained  by  adding  the  nuclear  repulsion  energy 

y 

and  reference  energy  to  the  approximate  AE  . In  the  absence  of 
analytic  methods,  the  derivative  of  the  total  energy  at  X=0  is 
computed  by  finite  difference  of  the  total  energy  evaluated  at  a small 
positive  and  a small  negative  value  of  X»  In  the  next  section,  we 
will  obtain  analytical  expressions  for  the  evaluation  of  ( 3AE/3X) | ^=q> 
so  that  together  with  analytic  derivatives  of  the  reference  energy  and 
the  nuclear  repulsion  energy,  we  will  have  the  analytic  derivative  of 
the  total  energy. 
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Figure  1.  The  Diagrammatic  Representation  of  H^CX). 


t°  {a  f i } = y 

t?f{atibtj}=  \f\J 

ttc{Qt|bVk}=  y y y 


Figure  2.  The  Diagrammatic  Representation  of  Coupled-Cluster 
Excitation  Operators  T^(X)>  ar>d  ^(X)* 
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Figure  3.  The  Diagrammatic  Representation  of  the  Coupled-Cluster 
Energy. 
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Figure  4.  The  Diagrammatic  Representation  of  the  T..(X)  Amplitude 
Equation  in  the  CCSDT  Model. 
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Figure  5.  The  Diagrammatic  Representation  of  the  T„(X)  Amplitude 
Equation  in  the  CCSDT  Model. 
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Figure  6.  The  Diagrammatic  Representation  of  the  T„(X)  Amplitude 
Equation  in  the  CCSDT  Model. 
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Figure  6 — continued 


FIRST  DERIVATIVES  OF  THE  COUPLED-CLUSTER  ENERGY 
Derivation  of  General  Expressions 

Let  us  return  to  Equation  (6),  where  the  Hamiltonian  is  in  normal- 
ordered  form  and  the  wavefunction  is  the  exponential  of  the  Coupled- 
Cluster  model. 


"T(X)  Hn(X)  eT(X)|0>  = AE(X)|0> 


(10) 


Let  us  define  the  following: 


3T(X) 


T = T(X) 


X=0 


9X 


X=0 


= ^N+  = ^N^X^  + ^N^X^ 


X=0 


HX  - fX+  WX  = 
HN  N N 


9fN(X)  9Wm(X)1 


N' 


l 9X  9X  J 


x=o 


AE  = AE(X) 


AEX  = 


9AE(X) 


x=o 


9X 


X=0 


Let  the  symbol  |P>  represent  the  space  of  orthonormal  determinants  in 
which  the  amplitudes  of  T(X)  and  AE(X)  are  determined.  Let  the  symbol 
| $>  represent  the  excited  determinants  only. 


1 


|PXP| 


(ID 


28 


Projection  on  the  left  by  <P|  and  evaluation  at  X=0  and  yields  the 
unperturbed  Coupled-Cluster  expressions. 


<0 1 (Hn  eT)c|0>  = AE  (12) 

and 

<*l  (Hn  eT)c|0>  = 0 (13) 

Of  course,  the  unperturbed  CC  model  (X=0)  is  a special  case  of  the 
perturbation-dependent  CC  model  discussed  in  the  previous  section. 
Equations  (12)  and  (13)  are  Equations  (8)  and  (9)  evaluated  at  X=0, 
respectively. 

The  derivative  of  Equation  (10)  with  respect  to  X,  evaluated  at 
X=0,  is  given  by  the  following,  where  we  have  used  the  fact  that 
excitation  operators  commute  with  one  another: 


[<Hn  eT)c,  TX|  ♦ (H*  eT)c 


|0> 


AEX|0> 


(14) 


Or  more  simply, 


((Hk 


eT)  TX)  + 
c c 


<HN  eT)c 


|0> 


AEX|0> 


V 

Projection  on  the  left  by  <P | yields  an  expression  for  AE  and  a set 
of  linear  equations  which  determine  the  derivative  t-ampl itudes : 


<0|((Hn  eT)cTX)c|0>  + <0|(HX  eT)c|0>  = AEX  (15) 

and 

<*I((Hn  eT)cTX)c|0>  + <*|(HX  eT)c|0>  = 


0 . 


(16) 
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For  the  CCSDT  model,  the  Equation  (16)  represents  three  coupled, 
linear  equations  for  the  TX,  TX  and  TX  amplitudes.  If  the  unperturbed 
CCSDT  amplitudes  are  available,  the  derivative  amplitude  equations  can 

Y 

be  constructed  and  solved.  The  derivative  energy  AE  can  then  be 
evaluated  by  Equation  (15). 

However,  it  is  not  necessary  to  solve  Equation  (16)  for  the 

V 

derivative  t-amplitudes  in  order  to  evaluate  AE  . Returning  to 
Equation  (14),  we  insert  1 = |PXP|  into  the  terms  involving  T and 
left  project  by  <P|  to  obtain 

<0|-TX|PXP|(HNeT)c|0>  + <0|(HXeT)c|0>  + <0 1 (HNeT)c | PXP |TX| 0>  - AEX 
and 

<*|-TX|PXP|(HNeT)c|0>  + <*|(HXeT)c|0>  + <*|  (HNeT)c|PXP|TX|0>  = 0 . 

But  TX  is  an  excitation  operator,  hence  <0|TX|P>=0.  Furthermore,  T is 

T 

the  solution  to  the  unperturbed  equation  <*|(H^e  ) c | 0>=0 . Also, 
recall  that  <0| (HNeT)c|0>=AE.  Now  we  have 

<0 1 (HX  eT)c|0>  + <0|(Hn  eT)c|*X*|TX|0>  = AEX 

and 

<$| -TX| 0>AE  + <*|(HX  eT)c|0>  + <*|(HN  eT)c | #><$| TX| 0>  = 0 . 

y 

The  second  equation  is  solved  for  <$|T  |0>: 

<*|TX|0>  = - <*|  (Hn  eT)c-  AE  |#>  <$|(HX  eT)c|0>  , 
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and  we  then  substitute  for  <$|TX|0>  in  the  AEX  expression. 

<0|(H^eT)c|0>  - <0|(HNeT)c|*>  <*|  («NeT)c-  AE  | * > <$| (HXeT)c | 0>  = AEX 

Now  we  define  an  antisymmetrized,  perturbation-independent,  de- 
excitation operator  A.  Let  A=A^+A2+A2+...,  where 


*1 


I A {i+a} 


if  a 


*2 


1 1J  + + 

- Z X {i+aj+b} 

4 ab 

if  j 

a,  b 


A 


n 


1 

(n!)2 


Z 


i j . . . 

X {i+aj+b. 
ab ... 


i f J f • • • 
a , b , . . . 


•} 


For  every  T in  T we  will  have  a corresponding  A in  A.  That  is,  for 
n n 

the  CCD  model  we  have  the  truncation  A = A^  For  the  CCSD  model  we 
have  A = A^  + A^  and  for  the  CCSDT  model  we  have  A = 

Ve  require  that  the  X-amplitudes  satisfy  the  condition: 


<0|A|*>  = - <0 1 (Hn  eT)c|*>  <*|  <Hn  eT)c-  AE  | t> 

V 

By  defining  A in  this  manner,  we  obtain  an  expression  for  AE  that  is 

V 

not  dependent  upon  T . We  have 
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<0|(HX  eT)c|0>  + <0|A|*X*|(HX  eT)c|0>  = AEX  , 

which  can  be  evaluated  if  we  first  solve 

<0 | A| #><*|  (Hn  eT)c-  AE  |*>  + <0 | (Hn  eT)c|$>  = 0 . 

Analogous  steps  have  been  often  used  in  recent  derivative  literature 
to  reduce  the  number  of  linear  equations  to  be  solved  from  3N  to  just 

one.  The  method  is  sometimes  referred  to  as  the  Z-vector  method  of 

u j 22,34 
Handy. 

We  now  wish  to  find  equivalent  forms  of  the  above  expressions 
which  can  be  readily  subjected  to  the  usual  diagrammatic  rules. 
Recalling  that  1 = |PXP|  = |0><0|  + |#X*|,  the  derivative  energy  can 
be  written  as: 

AEX  = <0|(HX  eT)c|0>  + <0|A  (HX  eT)c|0>  , (17) 

and  the  A equations  can  be  written  as 


<0|A  (HN  eT)cl$>  - AE  <0|A|*>  + <0|(HN  eT)c|$>  = 0 . 

If  we  consider 

<0|(Hn  eT)c  A|*>  = <0|(Hn  eT)c|PXP|A|*> 

= <0|(Hn  eT)c|0X0|A|$>  + <0|(Hn  eT)c  | *X*|  A|  *> 

= AE  <0 1 A | $>  + <0|(Hn  eT)c|*X$|A|$>  , 
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then  it  is  clear  that  ve  may  write  the  A equations  in  terms  of  the 

commutator  [A, (H„e  ) 1.  Of  the  commutator,  of  course,  only  the 
N C 

connected  product  is  non-zero.  Finally,  the  A equations  become: 

<0 | ( A (HNeT)c)cl*>  + <0|(HNeT)c|*><*|A|*>  + <0 | (HNeT)c | *>  = 0 . (18) 

Equations  (17)  and  (18)  summarize  the  CC  energy  derivative.  In 
Equation  (17),  AEX  is  expressed  in  terms  of  unperturbed  t-amplitudes, 
the  derivative  integrals  of  and  the  A-amplitudes.  Equation  (18)  is 
the  coupled,  linear  set  of  equations  which  determine  the  A-amplitudes. 
The  equations  are  general  in  that  they  apply  to  CC  models  of  any 
choice  of  truncation  of  T (and  the  corresponding  choice  of  A and 
projection  space  |P>).  Also,  the  equations  apply  to  CC  models  where 
non-Hartree-Fock  orbitals  are  chosen  as  the  spin-orbital  basis,  as 
well  as  when  SCF  orbitals  are  chosen.  The  equations  apply  generally 
to  CC  models  for  open-shell  systems,  where  Unrestricted  Hartree-Fock 
(UHF)  orbitals  are  used,  or  other  dif ferent-orbitals-for-dif ferent- 
spin  cases. 

y 

The  diagrammatic  representation  of  the  H*  operator  and  the  A 
operators  are  given  in  Figures  (7)  and  (8),  respectively.  The  fully- 
contracted  products  which  contribute  to  the  CCSDT  derivative  energy 
expression  are  represented  diagrammatically  in  Figure  (9).  The 
algebraic  expression  for  the  CCSDT  derivative  energy  is  given  in  Table 
(1);  each  term  corresponds  to  a diagram  in  Figure  (9)  and  is  so 
labeled. 

In  the  CCSDT  model,  Equation  (18)  represents  three  coupled,  linear 
equations  which  determine  the  amplitudes  of  the  A^,  A^,  and  A^ 
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operators: 


<0 | ( A (HNeT)c)c|S>  + <0|(HNeT)c|S>  = 0 

<0 1 ( A (HNeT)c)c|D>  + <0|(HNeT)c|SXS|A|D>  + <0 1 (HNeT)c  | D>  - 0 

<0 1 ( A (HNeT)c)c|T>  + <0|(HNeT)c|SXS|A|T>  + <0|  (HNeT)c|DXD|  A|T>  = 0. 

The  above  equations  are  represented  diagrammatically  in  Figures 

(10-12).  The  algebraic  interpretations  of  the  diagrammatic  equations 

are  given  in  Tables  (2-4),  respectively.  Note  that  although  there  are 

disconnected  parts  of  the  A^  and  A^  amplitude  equations,  there  are  no 

closed  and  disconnected  parts;  the  energy  derivative  is  necessarily 

linked.  For  the  CCSD  model  and  restricted  Hartree-Fock  (RHF) 

orbitals,  the  A^  and  equations  are  of  course  equivalent  to  that 

27 

obtained  by  Scheiner  et  al. 

In  Tables  (3)  and  (4)  we  use  the  symbol  P(occ)  (or  P(vir))  to 
represent  an  operator  that  generates  a sum  over  the  permutations  of 
occupied  (or  virtual)  labels  contained  in  parentheses,  with  each 
permutation  carrying  a parity  factor.  Specifically,  the  symbol  P(pq) 
found  in  Table  (3)  has  the  meaning: 

P(pq)  = identity  - interchange  p«q 

The  symbol  P(p/qr)  found  in  Table  (4)  has  the  meaning: 


P(p/qr) 


identity  - interchange  p«q  - interchange  p»r  , 
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where  the  slash  indicates  that  the  labels  q and  r are  actually 
equivalent — permutations  which  interchange  q and  r are  not  generated. 
The  symbol  P(vir|occ)  is  the  product  of  the  permutations  of  P(vir)  and 
P(occ).  For  example,  the  symbol  P(ab|ij),  found  in  Table  (3), 
generates  four  terms: 

P(ab|ij)  = identity  - interchange  a*b  - interchange  i*j 

+ interchange  a*b  and  i*j 

Finally,  the  symbol  P(a/bc | i/jk) , found  in  Table  (4)  generates  nine 
terms: 

P(a/bc | i/jk)  = identity  - interchange  a»b  - interchange  a»c 

- interchange  i»j  - interchange  i*k 
+ interchange  a*b  and  i*j  + interchange  a«b  and  i*k 
+ interchange  a*c  and  i»j  + interchange  a*c  and  i»k  . 

Comparison  of  the  X-amplitude  equations  with  the  unperturbed  t- 
amplitude  equations  (i.e.,  Figures  (4-6),  evaluated  at  X=0)  reveals 
that  the  X-amplitude  equations  possess  several  diagrams  which  are 
"upside-down"  versions  of  diagrams  in  the  t-amplitude  equations.  The 
algebraic  interpretations  of  such  diagrams  are  related  by  complex- 
conjugation;  if  the  spin-orbital  basis  functions  are  real,  the 
expressions  are  equivalent.  For  the  linearized  CC  models,  there  is  a 
one-to-one  correspondence  between  the  diagrams  in  the  X-  and  t- 
amplitude  equations;  hence,  the  X-amplitudes  are  equivalent  to  the 
unperturbed  t-amplitudes.  Consider  the  LCCD  model,  where  T=T2  and  the 
quadratic  T2  contributions  to  T2  are  not  permitted.  Specifically,  the 
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T2  equation  to  be  solved  consists  only  of  diagrams  (1)  and  (2a  e)  in 
Figure  (5).  The  corresponding  A2  equation  will  be  restricted  to 
diagrams  (1)  and  (8-12)  in  Figure  (11);  each  diagram  matches  a diagram 
in  the  T2  equation.  The  A equations  are,  as  already  mentioned,  linear 
equations;  thus,  for  the  CCD  model,  the  A2  equation  must  differ  from 
the  T2  equation  in  that  the  quadratic  diagrams  (5a-d)  in  Figure  (5) 
are  absent.  Instead,  we  find  the  similar,  related  diagrams  (21-27)  in 
Figure  (11).  If  we  denote  the  quadratic  contributions  to  T2  by  the 
function  Q*|?(I,  txt)  the  related  contributions  to  A^  are  given  by 
Qjj( t,  IxA)  + Q*J(t,  Axl). 

For  models  including  T^  we  have  previously  identified  the  linear 
contributions  of  Tj  to  T2,  represented  by  diagrams  (3a, b)  in  Figure 
(5);  analogously  there  are  linear  contributions  of  A^  to  A 2 — diagrams 
(4,5)  in  Figure  (11).  For  models  including  T^,  we  have  previously 
identified  the  linear  contributions  of  T^  to  T2,  represented  by 
diagrams  (4a-c)  in  Figure  (5);  analogously  there  are  linear 
contributions  of  A^  to  A^- diagrams  (33-35).  Consider  the  constant 
part  of  T^  and  linear  contributions  of  T^,  T2  and  T^  to  T^, 
represented  by  diagrams  (1),  (3a-c),  (2a-c)  and  (4)  in  Figure  (4). 
Analogous  diagrams  are  present  in  the  A^  equation  where  A^,  A^  and  A^ 
take  the  roles  of  T^ , T2  and  T3,  respectively.  The  analogous  diagrams 
are  diagrams  (1),  (3-5),  (25-27),  and  (78)  in  Figure  (10).  Finally, 
consider  the  linear  contributions  of  T2  and  T^  to  T^,  represented  by 
diagrams  (la,b)  and  (3a-e)  in  Figure  (6).  (Note:  there  are  no  linear 
contributions  of  Tx  to  T3-)  Analogous  diagrams  are  present  in  the  A3 
equation:  diagrams  (4,5)  and  (8-11)  in  Figure  (12). 

The  equivalence  of  A and  T in  the  linearized  CC  models  can  be  used 
to  advantage  in  the  evaluation  of  the  finite-order  MBPT  derivatives, 
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since  the  low-order  iterates  of  A and  T are  linear  and  therefore 
equivalent.  Only  the  fourth-order  contributions  from  quadruples  Q(4) 
will  require  special  consideration.  It  is  apparent  that  once  the 
equations  for  and  for  the  CCSD  model  have  been  programmed,  the 
A1  and  Aj  equations  for  the  CCSDT-1  model  will  follow  readily.  The 
same  routines  developed  for  inclusion  of  linear  triples  in  the  CCSDT-1 
t-amplitude  equations  can  be  employed  for  the  X-equations.  The  CCSDT- 
1 triples  routines  naturally  take  advantage  of  the  fact  that  since  T^ 
amplitudes  are  linearly  dependent  upon  the  T£  amplitudes,  it  is  not 
necessary  to  store  the  T^  amplitudes  at  any  time  during  the  CC 
iterations.  Rather,  the  back-contributions  of  T^  to  T2  and  to  T^  are 
computed  as  needed. 

In  order  to  evaluate  the  CCSDT  energy  expression  given  in  Table 
(1),  the  remaining  perturbation-dependent  parts  of  AE  — the 

X X 

contributions  arising  from  the  derivative  integrals  in  the  f^  and 
operators— must  be  made  explicit.  Recall  that  the  integrals  involve 
the  one-electron  spin-orbital  basis  functions  { 4>( X) } - Since  the 
orbital  basis  functions  are  usually  linear  combinations  of  some 
primitive  set  of  atomic  orbital  (AO)  basis  functions  {'KX)}>  the 
orbital  basis  functions  are  in  general  dependent  upon  X in  two  ways: 
the  primitive  AO  functions  may  be  dependent  upon  X;  the  coefficients 
may  be  dependent  upon  X-  Derivative  integral  programs  provide  the 
derivatives  of  one  and  two-electron  integrals  over  the  primitive  basis 
functions.  Coupled-Perturbed  Hartree-Fock  (CPHF)  theory  ’ provides 
the  derivatives  of  the  SCF  orbital  coefficients.  In  the  following 
sections,  we  exclusively  consider  the  choice  of  SCF  orbitals  as  the 


spin-orbital  basis. 
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Table  1. 

The  CCSDT  Derivative  Energy  Expression. 

AEX  = 

1 

5 

tjJ  <ij||ab>X 

[1] 

+ 

1 

2 

t*  tj  <ij | | ab>X 

[2] 

+ 

ta  f . X 
l la 

[3] 

+ 

X*  £ X 
a ai 

[1-1] 

- 

X1  f..X  ta 
a Ji  j 

[1-6] 

+ 

X1  f X tb 
a ab  i 
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+ 

X*  <ja| |bi>x  tj 
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Xj  f.  X 

b ia  ji 
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1 
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Xj  <ai | 1 bc>X  tjj 

H-3] 
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1 

2 

Xj  <jk|  | ia>X  tba 
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xj  f.  X tb  ta 

o ia  i j 
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+ 

X^  <ai | | bc>X  tb  t? 
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Xj  <jk|  | ia>X  tb  ta 

[ 1-14] 
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1 

? 

.k  , ,VX  _cab 

Xc  <ij 1 1 ab>A  tklj 
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Xk  <ij  | |ab>X  t“  tb 
c J 11  ki  j 
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+ 

1 

2 

.k  ...  i , , v X .ca  . b 
Xc<ij||ab>  t..  tk 

[1-9] 

+ 

1 

2 

xk  /..m  . v X *ab  . c 
Xc  <ij  1 1 ab>  tR.  tj 
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xk  ...  i i UvX  *a  . b . c 
X <ii  ab>  t,  t . t . 
c J 11  k l j 

[1-15] 

+ 

1 

5 

Xjj  <ab | | ij>X 

[2-1] 

1 

2 

X*j  <ab| |ci>X  tj 

[2-7] 

+ 

1 

2 

Xab  <kal li^>X  'k 

[2-8] 
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Table  1 — continued 
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1 
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<ai | 1 bc>X  t^b  t£ 

[2-18] 

- 

1 

5 

X3* 

ad 

<ai | | bc>X  tb£  tj 

[2-20] 

+ 

1 

2 

xf1 

be 

<jk| | ia>X  tbJ  ta 

[2-23] 

+ 

1 

5 

x;1 

be 

<jk| | ia>X  tb£  ta 

[2-21] 

- 

X*1 

be 

<jk| | ia>X  tba  t£ 

[2-19] 

- 

1 

2 

xk-> 

ad 

<ai | 1 bc>X  tb  tj  t? 

[2-36] 

+ 

1 

2 

x1! 

cb 

<3 k | | ia>X  t?  tb  ta 

[2-35] 
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Table  1 — continued 


+ 

1 

4 

Xdc  <« 

, Ksx  .dca  b 
lab>  tm  t. 

[2-27] 

+ 

1 

5 

Xdc  <ij 

, KSX  .dca  b 
lab>  tu.  tk 

[2-29] 

+ 

1 

5 

,1k  ..  . 
Xdc  <1J 

, X .dab  c 

lab>  tm  t. 

[2-28] 

+ 

1 

2 

Xcd  <1J 

. , X .ca  db 
|ab>  tk.  tXj 

[2-14] 

- 

1 

4 

Xcd  <« 

. , X .ca  db 
|ab>  t..  tkl 

[2-13] 

+ 

1 

15 

Xcd 

i lv X .ab  cd 
lab>  tkl  tjj 

[2-15] 

- 

1 

4 

XcJ 

|ab>  tk. 

[2-12] 

- 

1 

2 

Xcd  <1J 

, .SX  .cd  a b 
lab>  tk.  tj_  t. 

[2-33] 

+ 

1 

8 

Xcd 

,ahvX  .cd  a b 
lab>  t..  t1  tk 

[2-31] 

- 

1 

2 

x^  <ij 

1 ab>X  tS  tj  t) 

[2-34] 

- 

Xcd  <« 

1 ab>X  t“  tj  tb 

[2-32] 

+ 

1 

8 

»£  <‘J 

,auSX  .ab  d c 
Iab>  tkl  t.  t. 

[2-30] 

+ 

1 

4 

« <‘J 

,aKvX  .d  .c  a b 
lab>  t.  t.  tk  tx 

[2-37] 

+ 

1 

4 

xabd  <ab'  l“>X  *kj 

[3-1] 

- 

1 

4 

Xacb  <ka^  *lk 

[3-2] 

- 

1 

12 

>ilk  , X .abc 
abc  ji  llkj 

[3-9] 

+ 

1 

12 

\ijk  f X tdcb 
adc  ab  lijk 

[3-8] 

+ 

1 

24 

Xabe  <ab'  *'*>*  ‘“k 

[3-10] 

+ 

1 

24 

<UM*J>X  »SS 

[3-11] 

+ 

1 

4 

>W  <jal|bi>X  liy 

[3-12] 

+ 

1 

4 

Xabe  <abH'li>X  <Jk  'i 

[3-13] 
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Table  1 — continued 


1 

+ 4 

xijm 

“oca 

<kl|  |ij>X  tb‘  ta 

13-14] 

1 

' 2 

x ikl 
adc 

<ja | | bi>X  td?  tb 

[3-15] 

1 

" 2 

» ilk 
aed 

<ja| | bi>X  tjb  td 

[3-16] 

1 

' 12 

bed 

f X .bed  a 
xa  jki  ll 

[3-30] 

1 

" 12 

Jki 

Abcd 

f X tbca  td 
la  jki  l 

[3-31] 

1 

+ 12 

\ jki 
ade 

<ai | | bc>  ti 

[3-26] 

1 

" 4 

Jki 

ade 

<ai | | bc>X  tdab  tj 

[3-24] 

1 

+ 12 

\ jki 
ade 

<ai  | | bc>X  tdb‘  tj 

[3-28] 

1 

_ 12 

x ilm 
Dcd 

<jk|  |ia>X  tb'd  ta 

[3-27] 

1 

■ 12 

x ilm 
Abcd 

<jk | | ia>X  ta 

[3-29] 

1 

+ 5 

. ilm 
Dcd 

<jk|  |ia>X  td 

[3-25] 

1 

+ 4 

xkij 

cdb 

f X tca  tdb 
xa  lkl  ji 

[3-3] 

1 

+ 8 

xkij 

aed 

. mi.  xX  „ed  .be 
<ai  1 1 bc>  tkl.ty 

[3-6] 

1 

" 2 

kij 

ade 

x .ml  xX  ^db  . ec 
<ai||bc>  tH  tjt 

[3-4] 

1 

' 8 

x iml 
cdb 

✓ MM-  -xX  .cd  . ba 

<Jfcllla>  'jk  'ml 

[3-7] 

1 

+ 2 

x ilm 
cdb 

✓ -i  m • xX  .ca  . db 
<Jkl lla>  *lj  'ink 

[3-5] 

1 

+ 4 

Ml 

aed 

✓ • i it  xX  .ed  . b . c 
<ai  | | bc>  tRi  t^.  t1 

[3-34] 

1 

+ 2 

\kij 

aed 

✓ • i it  xX  .eb  . d c 
<ai | |bc>  tfcl  t.  tj 

[3-32] 

1 

“ 2 

x ilm 
cdb 

✓ •i  it-  xX  . c d . b .a 
<jk||ia>  tXj  tfc  tm 

[3-33] 

1 

* 4 

x iml 
cbd 

v X .ca  . b .d 
<3k| | ia>  tml  tj  tk 

[3-35] 

1 

" 24 

xklm 

ede 

✓ • • m u-xX  .cde  . ab 
<ijl!ab>  tku  tBj 

[3-19] 
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Table 


— continued 


1 

48 

.klm 

cde 

<ij 

| ab>X 

cde 

kij 

*2 

[3-22] 

1 

25 

.klm 

cde 

<ij 

| ab>X 

cda 

rklm 

‘I? 

[3-18] 

1 

5 

.klm 

cde 

<ij 

| ab>X 

cda 

rkli 

$ 

[3-17] 

1 

8 

.klm 

cde 

<ij 

| ab>X 

cda 

kij 

*£ 

[3-20] 

1 

48 

.klm 

cde 

<ij 

| ab>X 

cab 

rklm 

*« 

[3-23] 

1 

8 

.klm 

cde 

<ij 

| ab>X 

..cab 

rkli 

tde 

tmj 

[3-21] 

1 

12 

.klm 

cde 

<ij 

|ab>X 

cde 

rkli 

ta 

m 

•5 

[3-41] 

1 

25 

.kml 

cde 

<ij 

|ab>X 

cde 

lkij 

ta 

m 

[3-43] 

1 

12 

.klm 

cde 

<ij 

| ab>X 

cda 

^klm 

te 

l 

*5 

[3-42] 

1 

5 

.klm 

cde 

<ij 

| ab>X 

cda 

Kli 

te 

J 

< 

[3-45] 

1 

25 

.klm 

ced 

<ij 

| ab>X 

cab 

rklm 

t? 

l 

4 

[3-44] 

1 

5 

. lmk 
aec 

<ij 

|ab>X 

4l 

*2 

4 

[3-36] 

1 

2 

.1km 

aec 

<ij 

| ab>X 

41 

*3 

4 

[3-37] 

1 

8 

.lmk 

aec 

<ij 

| ab>X 

4l 

*2 

tc 

i 

[3-40] 

1 

8 

.mkl 

aec 

<ij 

| ab>X 

td? 

1J 

tca 

mk 

‘S 

[3-38] 

1 

2 

.lmk 

ace 

<ij 

| ab>X 

t? 

ki 

[3-39] 

1 

5 

.1km 

aec 

<ij 

| ab>X 

41 

t?  t 
J 

k ' 

i 

[3-46] 

1 

5 

.lmk 

ace 

<ij 

| ab>X 

4 : 

1 

e 

3 

4 

[3-47] 

Note:  Summation  over  all  labels  is  implied. 
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Table  2.  The  Equation  for  the  CCSDT  Derivative  Model. 
0 = fia  + <j  i I I ba>  t j 

+ Xb  fba 

Xj  f . . 
a ij 


1 

+ 2 
1 

+ 2 


1 

+ 2 
1 

+ 2 


xj  <ib| | aj> 
Xb  fi 


ja  *j 


XJ  f.. 
a lb 

Xj  <ci 


<cj 

Xb 

Xa  ^ 
*b  <ki 

Xb  <Jk 


XJ  <ki 
a 

Xj  <ki 


Xc  <kJ 

Xa 

Ji 


| ba>  t 


j 


|ba>  tj 

|ka>  tj 

|kb>  t$ 
|ca> 

be 
jk 

be 
jk 


| ca>  t 
| bc>  t 


c .b 


|ba>  tj  tj 


|ba>  t,'"  t T 


F tb 

k J 
|eb>  tj  tb 


. f i 1 7 
ca  cb  j 


lki  f tb 
Da  £jk  j 


Xjj  <cd | | ba>  tb 

<jMiki>  tj 

<jc||ak>  t!j 
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Table  2 


+ 


+ 


♦ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


— continued 


<ic|  |bk>  tj 


1 

2 

1 

2 

1 

2 

1 

2 


1 

2 


1 

5 

1 

2 

1 

2 

1 

5 


1 

2 


1 

2 


ca 


< j c | | bk>  tb 

xi1  fM 

ba  bj 

<bcl  la:i> 

<ib||jk> 

Ji  f tbc 

Abc  rka  jk 


f . t. 


be 


ba  ic  jk 

Xji  f tbc 
Aba  rkc  jk 

X31  £„.  t'  tb 
ca  kb  k j 


x3b  <di 

xdb  <dk 

xjk  <di 

da 

X*1  <dk 
da 

Xbj  <kl 

XU <3k 

^a  <ki 

<3k 

Xl1  <dk 
dc 

x5i  «* 

xdi 

Xcb  <k3 


be 

jk 

be 

jk 


. v _bc 
beX  tjk 


ca>  t 


ca>  t 


bc>  t 


be 

jk 


, . .be 

la>  t., 

Jk 

i v *bc 
la>  t., 

Jk 

i v -be 
lc>  *jk 

lc>  tbP 

Jk 

ba>  tP  tb 
k J 

CbX  tj  tj 

cb>  *k  *j 

lax  ,«  tb 
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Table  2 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


— continued 

xjj  <ki||lb>  ta  tb 

<kj||lb>  tb 

5 ^bc  <UI lda>  ‘jkl 

1 <kl||da>  tb'd 

5 ^ba  <xi 1 1 cd>  1 jkl 

5 Xba  <kl I I cd>  ‘jkl 

2 X«  aillbak  t'd 

2 Xcd  <lJHba>  *U  ‘j 

1 xJd  <klllba>  'kl  ‘j 

5 X^yiMdaX  tad  tj 

2 Xca  <jl||db>  ,ad  tb 

Xcb  'kl  'j' 

Xca  <li"db>  ‘kl  'j 

x«  <lj||db>  t‘d  tb 

2 <«Hdb>  ‘n  ‘S 

j Xkd  <ji| |cd>  t£d  tj 

2 >S<«llcd>  '«  ‘j 

2 Xdc  <lk I I ba>  td  t£  tb 

2 x$i  <m  lcb>  4 *£  ‘j 

Xb‘  <lj||cb>  td  t‘  tb 


1 

2 
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Table  2 — continued 


+ 

1 

2 

Xlji 

Dca 

fkl 

Jk 

+ 

1 

3 

xjki 

aeb 

<de 

ca> 

be 

*jk 

+ 

1 

8 

xjki 

aea 

<de 

bc> 

be 

jk 

+ 

1 

5 

bca 

<ki 

lm> 

be 

*jk 

+ 

1 

8 

xi" 

bca 

<jk 

lm> 

be 

*jk 

- 

1 

2 

Xlji 

Adbc 

<kd 

al> 

be 

ljk 

- 

1 

2 

xijk 

aba 

<id 

cl> 

be 

jk 

+ 

Xljl 

Adba 

<kd 

cl> 

+ 

1 

5 

\kji 

aea 

<de 

cb> 

+ 

1 

5 

x lmi 
cba 

<kj 

lm> 

*£ 

- 

Xlji 

Adca 

<kd 

bl> 

*£  *i 

- 

1 

2 

\kji 

cda 

<cd 

bk> 

<5 

3 

+ 

1 

2 

xkli 

cba 

<jc 

kl> 

tb 

J 

+ 

1 

5 

xjki 

bca 

<bc 

jk> 

+ 

1 

12 

xjki 

ebc 

<ei 

da> 

bed 

jkl 

- 

1 

5 

Jki 

ebc 

<el 

da> 

bed 
1 jkl 

+ 

1 

12 

xjki 

eba 

<ei 

cd> 

bed 

Jkl 

- 

1 

5 

xjki 

eba 

<el 

cd> 

bed 

Jkl 

- 

1 

12 

xmjk 

Abcd 

<li 

ma> 

bed 

Jkl 

- 

1 

12 

bed 

<kl 

ma> 

.bed 

Jkl 

+ 

1 

5 

xmjk 

bca 

<li 

md> 

bed 

*jkl 
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Table  2 — continued 


1 

12 

1 

12 

1 

3 

1 

2 

1 

2 

1 

2 

1 

2 

1 

2 

1 

2 


1 

3 

1 

2 

1 

2 

1 

2 

1 

3 

1 

2 


1 

2 

1 

2 

1 

12 


*bca  <kl""d> 

xjki  r bed 

^cd  rla  jkl 

xjkl  f bed 
Dca  rid  rjkl 

\jki  f tbcd 
“oca  rld  jkl 

xkji  f cd  b 
eda  rlb  rkl 

xkii  f cd  b 
eba  rjd  rkl  lj 


^jH  <el 

ecd 

,kli  y . 
Xecb  <eJ 

Xk«  <ei 
eca 

X <ej 
eca  J 

Xkjl  <el 
eca 

xkli  y . 

Xeba  <e^ 

xmki  . 
Xcdb  <1] 

X”kj  <11 
eda 

Xcda  <lj 


,i"Jl 

eda 


<kl 


x^a  <J1 

XSa  <« 

Xkr  <el 
eda 

X5ji  <lk 
dca 

|klj 

ede 


, . .cd  .b 
ba>  'ki  'j 

da>  t'd  tb 
db>  XU  ‘j 

db>  lii 

db>  tb 

„cd  .b 
cdX  tu  tj 

S «-ccl 

"a>  'kl  'j 

, v ..cd  .b 
"b>  'ki  'j 

mb>  tb 

kl  J 

mbX  tkl  tj 
mdX  tR1  t. 
mdX  'j 

cbX  tdx  1 1 t 

mbx  td  tk  t 
. . .cde  .b 

ba>  hfclm 


b 

j 

b 

j 
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Table  2 — continued 

1 

+ 12 

Xcde  <m^lba>  *klm  ') 

1 

+ 12 

X«*  <lm  1 1 ba>  tj 

1 

" 12 

»g;  <1*1 1->  'it: ') 

1 

+ 12 

Xcda  *klm  t5 

1 

~ 5 

Xcdb  <"lllea>  'klm  '] 

1 

' 5 

^ 'it: 

1 

+ ? 

Xcdi  <mj'leb>  'kS  'j 

1 

+ 5 

X«*  <lm|  [eb>  t£da 

1 

+ 12 

kcba  <jii|de>  *klm  *J 

1 

+ 5 

<"jnde>  'uS '] 

1 

+ 5 

&k  't:  *jk 

1 

‘ 2 

X«j  <mk  | | ca>  td* 

1 

+ 8 

>i£  <"1i'bc>  'S  'JE 

1 

" 5 

<«k||bc>  tda  tfk 

1 

+ 8 

Xjeb  <lml  lca>  'lm  'jk 

1 

+ 16 

Xj“  <lm|  |bc>  tda  tb‘ 

1 

" 2 

XJbl  <ki | | ec>  tda 

1 

' 5 

XX1  <jk|  | ec>  td*  tb‘ 

1 

+ 2 

xaii  <"kiiac>  'l:  ‘Je 

1 

+ 2 

X«J  <mk  | | ba>  ,da 

1 

+ 5 

Xdea  <mll lcb>  tda 

Table  2 


+ 


+ 


+ 


+ 


+ 


— continued 


1 

2 

> lki 
aea 

<mj  1 

| cb> 

*£ 

*£ 

tb 

3 

1 

8 

Jtji 

aea 

<lm| 

| cb> 

*£ 

tb 

J 

1 

3 

x lmi 
acb 

<kj| 

| ea> 

«£ 

*! 

1 

2 

\lmj 

Adca 

<ki  | 

| eb> 

lm 

«£ 

1 

2 

» lmi 
aca 

<kj| 

| eb> 

lm 

'K 

,b 

3 

Xlji 

aca 

<mk| 

| eb> 

de 

rlm 

*£ 

*s 

1 

8 

. lmi 
cba 

<kj  | 

| de> 

de 

rlm 

*£ 

1 

3 

xkji 

eda 

<ml  | 

| cb> 

te 

m 

'£ 

Note:  Summation  over  all  but  the  target  labels  i and 

implied. 
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Table  3.  The  Aj  Equation  for  the  CCSDT  Derivative  Model 
0 = <i j I I ab>  + P(ij |ab)  A*  f^b 

+ P(ij |ab)  A*  <kj | | cb>  t£ 

+ P(ij)  X^<cj||ab> 


P(ab) 

Ag  <i j 1 |kb> 

P(ij ) 

A^  <kj | | ab>  tj 

P(ab) 

Ak  <ij | |cb>  tj 

+ P(ij) 

Ak*  f.. 
ab  jk 

P(ab) 

Aij  « k 
ca  cb 

1 

+ 2 
1 

+ 2 

A* j <cd | | ab> 
<lj||ki> 

+ P(ab | ij ) 

<jc||bk> 

+ P(ab) 

Aij  f tc 
ca  rkb  rk 

+ P(ij) 

Aki  f tc 
ab  rjc  k 

xjj  <dk| | ab>  t£ 

+ P(ab | ij ) 

<djiicb>  «= 

+ P(ab) 

<dk||cb>  «' 

P(ab | ij ) 

x“  <kj||lb>  tj 

+ 

Xab  <‘5Hlc>  <k 

P(ij) 

X“  <kj||lc>  t' 

- j p(ij> 

Xcd  <«H«b>  'kl 

1 

+ 5 

X‘j  <kl  | | ab>  tj* 
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Table  3 — continued 


- 

1 

2 

P(ab) 

Xkl  <t]| 
ca  1 

db>  <bd 

+ 

P(ab | ij ) 

Xki  <lj| 
ca  J 1 

db>  lkdi 

- 

1 

2 

P(ab) 

Xij  <kl| 
ca  1 

db>  t'd 

+ 

1 

? 

Xab  <«  1 

cd>  tdd 

- 

1 

2 

P(ij) 

Xab  <«  1 

cd>  t'd 

+ 

1 

2 

X*j  <lk| 
dc  1 

ab>  td  td 

- 

P(ab | ij ) 

x$>  <u  i 

cb>  td  tj 

- 

P(ab) 

<ik' 

ebb  td  tj 

+ 

1 

2 

Xab  <4i 1 

dc>  'l  *k 

- 

P(ij ) 

X\l  <kJ 1 

dc>  td  t' 

+ 

Xkij  f 
cab  rck 

+ 

1 

2 

P(ab) 

Xcdl  <cd'|bk> 

- 

1 

2 

P(ij ) 

Xcab  <^Hkl> 

- 

cab  rkl 

*£ 

+ 

Xkij  f 
aab  rdc 

*£ 

- 

1 

2 

P(ab) 

Xdea  <de' 'cb>  'k 

- 

1 

2 

P(ij ) 

Xcab  *£ 

+ 

P(ab) 

xJcl  <kdi'bl>  i 

+ 

P(ij ) 

Xdab  <jd||cl>  'k 

+ 

<kd"ci>  'k 

+ 

1 

2 

P(ab) 

Xkij  f 
cda  rlb 

cd 

rkl 
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Table  3 — continued 


+ 

1 

2 

P(ij) 

.kli 

cab 

fjd 

cd 

ckl 

+ 

xkij 

cab 

fld 

cd 

ckl 

+ 

1 

2 

kij 

ecd 

<el 

ab> 

cd 

rkl 

- 

1 

2 

P(ab|ij) 

xk!i 

eca 

<ej 

| db> 

cd 

lkl 

- 

P(ab) 

xkij 

eca 

<el 

| db> 

cd 

rkl 

+ 

1 

5 

P(ij  ) 

»kli 

eab 

<ej 

| cd> 

cd 

rkl 

+ 

1 

2 

xkij 

eab 

<el 

| cd> 

cd 

rkl 

+ 

1 

2 

P(ab | i j ) 

»mki 

cda 

<lj 

| mb> 

cd 

rkl 

- 

1 

4 

P(ab) 

xmij 

cda 

<kl 

|mb> 

cd 

rkl 

- 

1 

2 

xtnkl 

cab 

<ij 

|md> 

cd 

rkl 

+ 

P(ij) 

. mki 
cab 

<lj 

| md> 

cd 

rkl 

- 

1 

2 

\mi  j 

cab 

<kl 

| md> 

cd 

rkl 

+ 

P(ab) 

xkij 

eda 

<el 

|.cb> 

d c 

h lk 

+ 

1 

2 

P(ij) 

xlki 

eab 

<ej 

| dc> 

d c 

*i  h 

+ 

Xlij 

eab 

<ek 

| dc> 

td  tc 
ri  \ 

- 

1 

2 

P(ab) 

x;ij 

dca 

<lk 

|mb> 

td  tc 
ll  lk 

- 

P(i  j ) 

xinki 

aab 

<lj 

| mc> 

td  tc 
rl  lk 

- 

xri 

dab 

<lk 

| mc> 

td  tc 
Z1  zk 

- 

1 

12 

P(ij) 

xkli 

cde 

<mj 

| ab> 

cde 

lklm 

+ 

1 

12 

\kij 

cde 

<lm 

| ab> 

cde 

Kim 

- 

1 

12 

P(ab) 

xklm 

cda 

<ij 

| eb> 

cde 

lklm 
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Table  3 — continued 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


5 P(ab|ij) 

XcJi  <*JMeb>  t'da 

\ P(ab) 

Xcdi  <it 

1 

12 

»SS <ij||de>  •£ 

l p(U) 

C 'it 

1 

5 

x^b  <i-i  i«*«>  xit 

1 

2 

Xdec  <mk||ab>  '£  *k 

2 P(ab | ij ) 

Xdea  ‘ill  'k 

2 P(ab) 

XJU  <mk|  | cb>  td* 

5 P(ab) 

Xdal  <1"l|cb>  'l»  'k 

2 P(ab | i j ) 

Xd"a  <ki"-b>  <t  'k 

1 

2 

Cb  <«!!«> 't  lk 

2 p(ij) 

Xdab  <«!!->  «£  'k 

P(ab) 

xJ‘i <mk'  teb>  ‘S  < 

P(ij) 

Xdab  <*"«>  tda  tj 
xdaj  <mk 1 1 ec>  lk 

1 

2 

Xdab  <1"H->  t' 

\ P(iJ) 

Xcab  <kjl|de>  'lm  'k 

1 

2 

X«b  <mk'|de>  '5m  *k 

2 P(ab) 

Xeda  <nll|cb>  *1  *1  'k 

2 P(U) 

Xeab  <"^|dc>  *m  'l  'k 
Xeab  <mk||dc>  t *1  'k 
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Note:  Summation  over  all  but  the  target  labels  i,  j,  a and  b is 

implied. 

P(pq)  is  an  operator  that  generates  two  terms: 

P(pq)  = identity  - interchange  p«q. 

P(ab|ij)  is  an  operator  that  generates  four  terms: 

P(ab|ij)  = identity  - interchange  a*b  - interchange  i*j 
+ interchange  both  a»b  and  i«j. 
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Table  4.  The  Equation  for  the  CCSDT  Derivative  Model. 


P(c/ab|k/ij) 

<ij | | ab>  Ak 

+ 

P(c/ab|k/ij) 

xil  <lkl ldc>  4 

+ 

P(c/ab|k/ij ) 

Xab  fkc 

- 

P(a/bc|k/ij) 

<dk| |bc> 

+ 

P(c/ab | i/jk) 

Xab  llc> 

+ 

P(a/bc|k/ij ) 

<lk| |bc>  td 

+ 

P(c/ab | i/jk) 

<jk||dc>  tj 

- 

P(k/ij) 

Xlij  f 
abc  rkl 

+ 

P(c/ab) 

Xijk  f 
Adab  dc 

+ 

1 

2 

P(a/bc) 

Xdea  <deMbc> 

+ 

1 

2 

P(i/jk) 

Xabc  <jkMlm> 

+ 

P(c/ab|k/ij) 

X]H  <kd||cl> 

- 

P(c/ab) 

Xijk  f td 
Adab  rlc  ll 

- 

P(k/ij) 

Xlij  f td 
abc  rkd 

- 

P(a/bc) 

Xeda  <el"bc>  4 

+ 

P(c/ab | k/ij ) 

Xeab  <ek| |dc>  *1 

- 

P(c/ab) 

Xeib  <ell'dc>  4 

- 

P(c/ab|k/ij ) 

Xdab  <lk| |mc>  4 

+ 

P(i/jk) 

Xabc  tmd>  4 

+ 

P(k/ij) 

Xabc  <lkHn,d>  4 

+ 

1 

2 

P(a/bc|k/ij) 

Xdea  <"kHbc>  'lm 
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Table  4 — continued 

* J PC'bc)  Xjij  <l"Hbc>  t£ 

+ \ P(c/ab|i/jk)  xj“  <jk||ec>  tj' 

* P(c/ab|k/ij ) X^j  <mk||ec>  tj* 

+ j P(c/ab)  Xdib  <lm||ec> 

* J PCi/3k)  X“  <jk|  |da> 

* \ P(k/ij)  <mk!  |de>  tj* 

* l P(a/bc)  X*£  <ml||bc>  t'  t$ 

P(c/ab|k/ij)  X“J  <mk||dc>  t'  tj 

+ P(c/ab)  X^  <ml||dc>  t'  tj 

* 2 PU/JX)  xj“  <jk||ed>  t*  t* 

* PCk/ij)  <lk| |ad>  t*  t{ 

Note:  Summation  over  all  but  the  target  labels  i,  j » k,  a,  b and  c is 

implied. 

P(p/qr)  is  an  operator  that  generates  three  terms: 

P(p/qr)  = identity  - interchange  p«q  - interchange  p»r. 

P(a/bc | i/jk)  is  an  operator  that  generates  nine  terms: 

P(a/bc  | i/jk)  = identity  - interchange  a«b  - interchange  a*c 
- interchange  i*j  - interchange  i*k 
+ interchange  both  a*b  and  i«j  + interchange  both  a«b  and  i»k 

+ interchange  both  a»c  and  i*j  + interchange  both  a*c  and  i»k. 
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Expressions  for  the  Choice  of  SCF  Reference 

For  canonical  SCF  orbitals,  all  terms  in  Table  (1)  involving  the 
derivative  integrals  of  fX  vanish,  except  for  those  involving  the 
diagonal  elements  f..X  = e.X  and  f X = e X— the  derivative  SCF 
eigenvalues.  The  results  of  CPHF  theory15  are  now  summarized  for 
canonical  SCF  orbitals  { 4*( X) } expanded  in  a primitive  AO  basis  set 
{*(X)}. 

In  this  section,  all  orbitals  are  assumed  to  be  unperturbed;  i.e. 
evaluated  at  X=0,  unless  otherwise  shown  as  a function  of  X- 


Derivative 

orbitals, 

evaluated  at  X=0,  are 

indicated 

by  the  super- 

script  X- 

I p(X)> 

= ♦pCX) 

Ip>  = Ip(x)>Ix=0 

Ipx> 

3Ip(X)> 

ax 

|y(x)> 

■ v*> 

|y>  = 1 u(x)> lx=0 

|yx> 

3|y(x)> 

ax 

As  usual,  the  labels  a,b,...  and  i,j,...  are  used  to  represent  the  SCF 
virtual  orbitals  and  occupied  orbitals,  respectively,  while  p,q,r,... 
are  unrestricted.  The  labels  y,v, X, a, ...  are  used  to  represent  the 
primitive  AO  basis  functions.  The  derivative  of  an  SCF  orbital, 
evaluated  at  X=0,  is  given  by  a linear  combination  of  unperturbed  SCF 
orbitals  and  a term  arising  from  the  derivative  of  the  primitive  AO 
basis  functions: 

| aX>  = l Ufa|f>  + £ \alk>  + 1 Cpa|yX>  (19) 

|iX>  = E UX|f>  + I UX|k>  + ZC.|yX>  . 

f 1 k*i  w 


(20) 
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The  SCF  coefficients  for  the  unperturbed  molecular  orbital  p are  given 
bv  C . The  coefficients  UX  are  called  the  CPHF  coefficients.  The 

3 yp  pq 

virtual— occupied  block  of  the  CPHF  coefficients  are  determined  by 
solving  the  linear  equation: 


E (A' 
bj 


ai , bj 


(21) 


where  A i <ab I I i J > + 1 | ib> 

ai , bj  = e - t. 


and  e is  the  unperturbed  orbital  eigenvalue  of  orbital  p.  Once  the 
P 

virtual-occupied  block  is  determined,  the  occupied-occupied  and 
virtual-virtual  blocks  are  then  determined  as  well: 


,x 

-1 

QX  + E UX  (<km 

kl  „ m 
g>m 

I ig>  + <kg 1 

|im>)  ] 

(22) 

'ki  = 

(Vci> 

lx 

-1 

[ QX  + E UX  (<fm| 

fa  n-  m g"1 

g»m 

1 ag>  + <fg| 

| am> ) ] 

(23) 

'fa  = 

<crEa> 

The  occupied-virtual  CPHF  coefficients  are  determined  by  the 
derivative  of  the  SCF  orthonormality  condition: 


V * 

(U*  ) + 

qp 


o 


(24) 


The  elements  of  the  general  one-electron  perturbation  matrix  Q 


X 


are 


given  by: 


SX  t 

pq  pq 


Z SX  <pl | | qk> 
k,l 


+ 


Z 

y,v 
X,  a 


(V*% 


9<yX|  | va> 

ACT  3X 


X=0 


where 


2 (8 


+ e ) 

p q 


when  p and  q are  both 
virtual  or  both  occupied 


(25) 


v. 


otherwise 
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The  symbol  P»  represents  an  element  of  the  SCF  density  matrix  in  the 
AO  basis.  The  matrix  hXq  consists  of  the  derivative  of  the  matrix 
elements  of  the  one— electron  part  of  the  Hamiltonian  in  the  primitive, 
perturbation-dependent  basis  {<KX)b  but  transformed  into  the 
unperturbed  orbital  basis  {4>). 


E 

y,  v 


a<u(x)lh(x)lv(x)> 

3X 


x=0 


Similarly,  the  matrix  SXq  is  the  derivative  of  the  overlap  matrix  of 
the  perturbation-dependent  primitive  basis  functions,  transformed  into 


the  unperturbed  orbital  basis. 


E 

y,v 


(Cyp^  Cvq 


3<y(X)|v(X» 

3X 


X=0 


(26) 


For  a uniform  external  field  perturbation  in  the  z-direction,  the 


elements  of  the  perturbation  matrix  Q are  simply: 


For  the  delta-function  perturbation 


X 

pq 

36 


= <p|z|q>  • 

located  at  the  point  r,  ve  have: 


QX  = hX  = <p| S(r'-r) |q>  = •Mr)  4>(r)  • (27) 

pq  pq  p q 

The  derivative  of  the  SCF  orbital  eigenvalues,  the  diagonal 
elements  of  the  Fock  matrix,  are  given  by: 


eX  = QX  + E UX  (<am||ag>  + <ag||am>) 
a aa  gm 

eX  = q£.  + E UXm(<im | | ig>  + <ig||im>)  . 

1 1 g,m  g 

Now  ve  could  use  the  above  identities  to  develop  the  expression 
for  AEX,  except  that  the  occupied-occupied  and  virtual-virtual  CPHF 
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coefficients  (Equations  (22)  and  (23))  possess  denominators  involving 
the  difference  of  unperturbed  occupied  or  virtual  eigenvalues.  Such 
denominators  must  be  eliminated  by  subsequent  factorization  of  the 
derivative  energy  expression;  otherwise  the  expression  will  not  be 
applicable  to  molecular  systems  with  degeneracies  or  near 
degeneracies.  Since  CC  theory  and  its  finite-order  MBPT 
approximations  are  invariant  under  unitary  transformations  so  long  as 
occupied  and  virtual  orbitals  are  not  mixed,  we  may  make  a particular 
choice  of  non-canonical  SCF  orbitals  to  simplify  the  derivative 
expressions  without  loss  of  generality.  Rather  than  introduce 
singularies  which  must  be  later  eliminated,  we  follow  the  development 
of  Handy  et  al.  for  MBPT(2)  gradients29  by  choosing  non-canonical  SCF 
orbitals  such  that: 


U£a  ’ ' 2 

Sfa'£> 

(28) 

“ki  ' - 2 

Ski  1 k>  • 

(29) 

. 37 

as  found  in  early  studies  by  Moccia 

and  Pulay.3-^ 

The  problematic 

denominators  of  the  occupied-occupied  and  virtual-virtual  CPHF 

coefficients  for  canonical  orbitals  are  not  present;  the  virtual- 

occupied  CPHF  coefficients  are  determined  as  in  Equation  (21).  For 

the  choice  of  non-canonical  orbitals,  however,  the  off-diagonal 

elements  f ,X  and  f . .X  are  not  zero.  These  matrix  elements  are  given 
ab  ij 

by: 


f X = QX  + E UX  (<am  bg>  + <ag  |bm>)  (30) 

ab  ab  „ „ gm 

g>m 


f..X=  QX.  + Z UX (<im||jg>  + <ig||jm>)  . 

1J  13  g,m  g 


(31) 
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Equations  (28-31)  and  Equation  (24)  are  to  be  applied  to  the  CCSDT 
energy  derivative  expression  given  in  Table  (1).  First,  gathering 
together  the  common  derivative  two-electron  integrals  and  the  non-zero 
derivative  Fock  matrix  elements,  we  define  the  intermediates  I"  and  D 
for  a more  compact  expression: 


AEX 


+ 


T(ij ,ab)  <ij | | ab>X 

r(kl,ij)  <kl | | ij>X 

T(ai,bc)  <ai||bc>X 


r(ab,cd)  <ab||cd>X 
r( ja, bi ) <ja | |bi>X 
T(jk,ia)  <jk | | ia>X 


where  summation  over  all  labels  is  implied.  The  D and  T intermediates 
for  the  CCSDT  model  are  defined  in  Tables  (5)  and  (6),  respectively. 
For  the  CCSDT-1  model,  the  contributions  to  the  T's  which  arise  from 
terms  [2-27],  [2-28]  and  [2-29]  are  excluded.  In  addition,  the 
contributions  from  terms  involving  are  excluded,  with  the  exception 
of  terms  [3-1],  [3-2],  [3-8]  and  [3-9].  We  will  show  that  the  D 
intermediates  so  defined  are  the  occupied-occupied  and  virtual-virtual 
blocks  of  the  density  matrix. 

Now,  using  the  identities  from  CPHF  theory,  we  expand  the 
derivative  integrals,  and  obtain  the  following: 


&EX  = D..  (QX.  + UX  (<ie||jm>  + <im| |je>)) 

lj  ij  em 

+ D , (QX,  + UX  (<ae||bm>  + <am||be>)) 
ab  ab  em 

+ <r(ij,ab)  - r(ji,ab))  (UX.<ej||ab>  - 1/2  SX<mj||ab>) 


(32) 
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+ (T(ij ,ab)  - T(ij , ba) ) C [ -U S ] <i j | |mb>  - 1/2  S*a<ij||eb>) 

+ (T(ab,cd)  - T(ba,cd)  + T(cd,ab)  - I\cd,ba)) 

x ^-Uam-S^m]<mb||cd>  ~ 1/2  Sea<eb||cd>) 

+ (T(ki,ij)  - r(ik,ij)  + r(ij,ki)  - r(ij,ik)) 
x (Ugjc<el  | | ij>  - 1/2  sjk<ml||ij>) 

+ (T(ja,bi)  + T(ib,aj))  (U^<ea||bi>  - 1/2  S^<ma||bi>) 

+ (I\ja,bi)  + r(ib,aj))  U-Uam~Sam]<jm|  |bi>  ' 1/2  Sea<je||bi>) 
+ IXai.bc)  U-Ugm-S^m]<mi|  |bc>  - 1/2  S*a<ei||bc>) 

+ T(ai , be)  (U^.<ae||bc>  - 1/2  S*.<am||bc>) 

+ (T(ai ,bc)  - T(ai , cb) ) ( [-U^-S^]  <ai  | | mc>  - 1/2  S*b<ai||ec>) 
+ (I\ jk, ia)  - T(kj,ia))  (U^<ek||ia>  -1/2  S^<mk||ia>) 

+ r(jk,ia)  (Ugi<jk|  | ea>  - 1/2  S*.<jk|  |ma>) 

+ r(jk,ia)  ( [-U^m-S^m]<jk| |im>  - 1/2  S*a<jk||ie>) 

+ AO  derivative  part 

Contributions  arising  from  the  derivative  of  the  primitive  AO  basis 

functions  when  Equations  (19)  and  (20)  are  applied  to  the  derivative 

two-electron  integrals  are  called  the  "AO  derivative  part".  The  AO 

derivative  part  is  not  explicitly  shown  here  because  the  following 

discussion  is  concerned  with  manipulations  of  the  other  terms. 

Rather,  the  AO  derivative  part  is  shown  in  final,  computational  form 

26 

in  Table  (7).  In  Table  (7),  we  follow  Fitzgerald  et  al.  and  define 
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a perturbation-independent  quantity  r(yvptr)  by  the  application  of  MO- 
to-AO  transformations  to  the  T intermediates.  The  transformation  of 
3N  sets  of  derivative  integrals  to  the  MO  basis  is  thereby  avoided. 
Instead,  only  the  one  transformation  in  the  construction  of  the 
perturbation-independent  r(yvpo)  is  necessary.  The  AO  derivative  part 
of  each  of  the  3N  gradients  is  computed  by  the  dot-product  of  r(yvpa) 
with  the  AO  derivative  integrals  as  they  are  generated  for  each  of  the 
3N  degrees  of  freedom. 

We  gather  together  the  common  terms  of  Equation  (32)  by  defining 
computational  intermediates  I and  X such  that: 


AE' 


‘ * 2 Dab  Qab  * 2 C.  Xai  U«i 

a,  i 


j 


1 3 iJ 


a,  b 


I I,.  S?.  + I I , SX,  + 2 E I.  S*. 


1J  13 


ab  ab 


a,  i 


X 

ai  ai 


i»J  “ a,b 

+ AO  derivative  part 

The  intermediates  X and  I are  summarized  in  Table  (8).  Now  recall 
that  the  CPHF  coefficients  are  the  solutions  to  a linear  equation. 
Hence,  we  use  the  Z-vector  method  ’ to  eliminate  the  CPHF 
coefficients  from  the  expression: 


Q 


2 1 X . UA. 
. ai  ai 
a,  i 


M. 


■ 2 E.  Xai  E.<A‘\i,bj  c.  -€h 

ai  bj  J b 


= 2 1 .Dbj  QSi  * 

b,j  J 


where  the  summation  over  the  labels  a and  i has  been  carried  out  in 
the  term  at  the  far  right.  The  virtual-occupied  block  of  the  matrix 
is  the  solution  to  the  perturbation-independent  linear  equation: 


l D..(e.  - e.)  A..  . = X . 

bj  bj v j b'  bj,ai  ai 


(33) 
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Finally  we  have  the  computational  form  of  the  energy  derivative: 


AEX  = ED..  QX. 

* E Dab°ab  * 

2 E.  Dai 
a,  i 

QX. 

ai 

i > j 1J  1J 

a,  b 

+ E I. . SX. 

. . n n 

* E Eab  Sab  * 

2 E.  Jai 
a,i 

SX 

ai 

i » j 3 3 

a,  b 

+ AO  derivative  part 

In  the  special  case  of  the  applied  external  field,  Equation  (34) 
reduces  to 


A EX 


E 

p.q 


Thus,  for  the  delta-function  perturbation  (see  Equation  (27))  we  have 


“x  * p1,  Dpq  *p(t)  Vr)  ’ 

and  we  identify  the  matrix  D as  the  matrix  representation  of  the 

correlation  correction  to  the  electron  density.  When  added  to  the  SCF 

density  matrix,  the  total  response  density  matrix  is  obtained.  Since 

the  molecular  orbitals  are  allowed  to  "relax"  in  the  presence  of  the 

perturbation,  the  response  density  necessarily  contains  all  orbital 

relaxation  effects,  as  opposed  to  a density  constructed  by  taking  the 

expectation  value  over  a given  unperturbed  wavefunction.  It  is  well 

known  that  iteration  of  the  CPHF  equation  (Equation  (21))  and 

subsequent  substitution  into  the  SCF  second  derivative  energy 

expression  yields  certain  classes  of  double-perturbation  diagrams  to 

all  orders. ^ Similarly  the  first  derivative  of  the  MBPT/CC  energy  is 

39 

actually  an  infinite  sum  of  diagrams  in  double  perturbation  theory. 
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Table  5.  The  Occupied-Occupied  and  Virtual-Virtual  Blocks  of  the 
CCSDT  Density  Matrix. 


D.  . = - 

ij 

Xj  ta 
a l 

1 

+ 2 

x»i 

ab  ki 

1 

" 12 

.ilk  .abc 
Xabc  'lki 

Dab  " 

X*  ,b 
a l 

1 

‘ 2 

X1J 

ac  ij 

1 

+ 12 

.ijk  deb 
adc  ljk 

Note:  Summation  over  all  but  target  labels  is  implied. 


78 


Table  6. 
T(ij ,ab) 


The  T Intermediates. 


1 

Xij 

3 

ab 

1 

. ab 

+ 

3 

t . . 
il 

1 

.a  .b 

+ 

2 

1 

xk  _cab 

+ 

3 

Xc  Vij 

.k  _ca  „b 

+ 

Xc  Vi  'j 

1 

,k  „ca  .b 

+ 

2 

Xc  *ij  V 

1 

Nk  ab  .c 

+ 

2 

X t,  . t . 
c ki  j 

>k  a „b  .c 

+ 

X t,  t.  t. 
c k l j 

1 

.Ik  dca  b 

+ 

3 

ac  rlki  j 

1 

xlk  dca  b 

+ 

3 

Adc  rlij  rk 

1 

.Ik  dab  c 

+ 

3 

ac  llki  j 

1 

■xkl  ca  db 

+ 

2 

cd  rki  rlj 

1 

xkl  ca  db 
cd  lj  rkl 

3 

1 

xkl  ab  cd 

+ 

16 

cd  rkl  lj 

1 

xkl  ab  cd 

3 

cd  tki  rlj 

1 

xki  cd  a 
cd  rki  rl 

2 

1 

xlk  _cd  a 

+ 

8 

Xcd  *ij  'l 

1 

xkl  ca  d 

2 

Xcd  lkl  li 

xkl  ca  d 
cd  *ki  j 

1 

xkl  .ab  .d 

+ 

8 

Xdc  lkl  li 

1 

xkl  d tc  . 

+ 

3 

X.  t . t . t 
dc  l j 

J 

i 

tb 

j 

tb 


t? 

3 

a 

k 


t 


b 

1 
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Table  6 — continued 


1 

.klm 

.cde 

ab 

25 

cde 

hli 

hj 

1 

.klm 

..cde 

ab 

+ 

48 

cde 

hij 

hm 

1 

.klm 

cda 

eb 

25 

cde 

lklm 

hj 

1 

. klm 

_cda 

eb 

+ 

5 

cde 

hli 

hj 

1 

.klm 

cda 

hij 

eb 

hm 

8 

cde 

1 

.klm 

cab 

de 

+ 

48 

cde 

lklm 

hj 

1 

» klm 

_cab 

de 

8 

cde 

tkli 

hj 

1 

.klm 

..cde 

ta 

tb 

3 

12 

cde 

hli 

h 

1 

.kml 

_cde 

ta 

tb 

+ 

25 

cde 

hij 

t 

m 

h 

1 

.klm 

cda 

fe 

tb 

12 

cde 

him 

h 

h 

1 

.klm 

..cda 

fe 

tb 

5 

cde 

tkli 

h 

h 

1 

.klm 

cab 

te 

td 

+ 

25 

ced 

tklm 

l . 
1 

J 

1 

x lmk 

de 

ca 

tb 

+ 

5 

aec 

hi 

mk 

h 

1 

.1km 

de 

ca 

tb 

2 

aec 

hi 

hj 

h 

1 

.lmk 

*s 

ab 

t? 

3 

+ 

8 

aec 

mk 

1 

.mkl 

de 

ca 

tb 

+ 

8 

aec 

t . . 
ij 

mk 

h 

1 

.lmk 

da 

cb 

te 

2 

ace 

hm 

hi 

J 

1 

.1km 

4! 

a 

+ 

5 

<iec 

L . I 

J 

k 

1 

.lmk 

da 

t?  t® 
i J 

+ 

5 

clce 

hm 
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Table  6- 


r(ab,cd) 


r(kl,ij) 


r(ja,bi) 


T(ai ,bc) 


continued 


1 

8 

1 

+ 5 
1 

+ 25 
1 

+ 5 


1 

8 

1 

+ 5 
1 

+ 25 
1 

+ 5 


1 

+ 5 
1 

" 2 

1 

‘ 2 


1 

2 


1 

' 5 
1 

" 2 


1 

5 


ab  13 

ab  j 1 

x ijk  ecd 
abe  ljk 

xjki  ec 
abe  rjk  r 


X!j  t?b 
ab  kl 


xbi  ik 


,ijm  abe 
abe  mkl 

X*jm  tb?  t 
bca  mk 


d 

i 


a 

1 


t 


b 

j 


xik  tc 

ac  3 


t 


b 

k 


ilk  deb 
adc  rlkj 

ikl  dc  b 
adc  lkj  ll 

ilk  cb  d 
acd  rlk  3 


tT 
a 3 


\kl  t 

ad 

XJ*  t 
ad 

x3i  t 

ad 

Xjk  t 
ad 


dbc 

kji 

db 

jk 

db 

ji 

be 

jk 


t 

t 

t 


c 

i 

c 

k 

d 

i 


81 


Table  6- 


T(jk,ia) 


■continued 


1 

2 

X«  tS  t?  t? 
ad  k l j 

1 

\ j kl  tdeb  c 

12 

ade  jkl  i 

1 

\ jkl  .deb  c 

4 

ade  rjki  Z1 

1 

\ jkl  tdbc  e 

12 

ade  ljkl 

1 

\klj  .ed  be 
aed  rki  rlj 

8 

1 

\klj  .db  ec 

2 

ade  rkl  ji 

1 

\kjl  ted  b c 

5 

aed  rki  rl 

1 

XklJ  t'?  t?  tc 
aed  kl  l j 

2 

1 

2 

be  j 

1 

5 

x*kj  t?- 

oed  kj 

1 

X* 

b jk 

2 

X*  t1?  t? 
b j k 

1 

.il  eba 

5 

cb  tljk 

1 

\il  tbc  a 

2 

V 'ij  'k 

1 

x?1  t? 

be  jk  1 

5 

xil  .ba  c 

X,  t,  . t. 
be  lj  k 

1 

. il  r „b  _a 

2 

X , t . t,  t, 
cb  j k 1 

1 

>ilm  bed  a 

12 

Abcd  rlmj  Tk 

1 

.ilm  bed  a 

12 

Abcd  tljk  m 

1 

xilm  bca  d 

5 

Abcd  rlmj  rk 

1 

.iml  cd  ba 

8 

cdb  jk  ml 
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Table  6 — continued 


,ilm  ca  db 
cdb  rlj  mk 

.ilm  cd  b a 
cdb  rlj  lk  m 

ximl  ca  b d 
cbd  ml  j :k 


& t* 
ab  1 


1 

S 


t 


cb 

li 


Note:  Summation  over  all  but  the  target  labels  is  implied. 
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Table  7.  The  AO  Derivative  Part  of  the  Coupled-Cluster  and  MBPT 
Energy  Derivatives. 


AO  derivative  part  = E <p(X) v(X)  | p(X) <KX)>*  r(pvpa)  , 

U,v, 
p,  a 


where 


rojvpc)  . I r<ij,ab)  (Cca  c3b  - ccb  cre)  C*J  C*J 
a,b 

* £ r(ab, cd)  (Cpc  c5d  - Cpd  Cpc)  c*a  <b 

a,  b 
c,d 


* * 


+ I r<ij,ki)  (Cpk  Cpl  - Cpl  cok)  cui  Cv. 

* y J 
k,l 

* t <c„b  Cci  - Cpi  Cpb>  Cbj  C 

i , a 

j, b 

+ E r(ai , be)  (C  Cffc  - C Cab)  C*a  <. 
i,a 
b,  c 

+ E T(  jk,  ia)  (C  Caa 
i,a  • 


* * 


C C . ) C . c , 
pa  01'  VJ  vk 


j»k 


Note:  <yv| pa>X  denotes  the  derivative  of  an  ordinary  AO  two-electron 

integral,  evaluated  at  X=0. 
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Table  8.  The  I and  X Intermediates  Used  in  the  Evaluation  of  CCSDT 
and  MBPT  Energy  Derivatives. 

Ii.  = -1/2  (T(jk,ab)  - T(kj,ab))  <ik||ab> 

-1/2  (T(jl,mk)  - r(lj,mk)  + T(mk,jl)  - r(mk,lj>)  <il| | mk> 
-1/2  (T(ja,bk)  + T(kb,aj))  <ia||bk> 

-1/2  T(aj ,bc)  <ai| |bc> 

-1/2  (T(jk,la)  - T(kj , la) ) <ik||la> 

-1/2  T(lk, ja)  <lk| | ia> 

Iab  = -1/2  (r(ij,bc)  - r(ij,cb))  <ij | | ac> 

-1/2  (T(be,cd)  - T(eb,cd)  + T(cd,be)  - T(cd,eb))  <ae||cd> 
-1/2  (T(jb,ci)  + T(ic,bj ))  <ja||ci> 

-1/2  T(bi,cd)  <ai| | cd> 

-1/2  (T(di ,bc)  - T(di , cb) ) <di||ac> 

-1/2  T(jk,ib)  <jk| | ia> 

I . = -1/2  (T(kj,ab)  - T(kj , ba) ) <kj||ib> 

ai 

-1/2  (T(ab, cd)  - T(ba,cd)  + T(cd,ab)  - T(cd,ba))  <ib||cd> 
-1/2  (T(ja,bk)  + T(kb,aj))  <ji||bk> 

-1/2  T(aj ,bc)  <ij | | bc> 

-1/2  (T(dj ,ac)  - T ( d j , ca) ) <dj||ic> 

-1/2  I\jk,la)  < j k | | li> 
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Table  8 — continued 


+1/2  Dkj  (<ka | | ji>  + <ki | |ja>) 

+1/2  Dcb  (<ca| | bi>  + <ci t | ba> ) 

+1/2  (r(ij,cb)  - I\ji,cb))  <aj||cb> 

+1/2  ( r ( i 1 , k.  j ) - r(li,kj)  + T(kj,il)  - T(kj,li))  <al||kj> 
+1/2  (T(ic,bj)  + r(jb,ci)>  <ac||bj> 

+1/2  T(di , be)  <da| |bc> 

+1/2  (T(ik,jb)  - T(ki,jb))  <ak| | jb> 

+1/2  T( jk, ib)  <j k | | ab> 

Note:  Summation  over  all  but  the  target  labels  is  implied. 

The  occupied-occupied  and  virtual-virtual  blocks  of  the  CCSDT 
density  matrix  D are  given  in  Table  (5). 

The  T intermediates  are  defined  in  Table  (6). 
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MBPT  Energy  Derivatives 

Now  we  extract  the  finite-order  MBPT  derivative  from  the  CCSDT 
derivative  energy  expression.  Inspection  of  Equation  (34)  reveals 
that  for  the  derivative  of  the  n^-order  correlation  energy,  the  n 
order  MBPT  density  matrix  D and  the  nth-order  I intermediates  are 
required.  We  specify  the  construction  of  the  occupied-occupied  and 
virtual-virtual  blocks  of  the  density  matrix  D and  the  T intermediates 
appropriate  for  various  MBPT  models  in  Tables  (9)  through  (12).  The  X 
and  I intermediates  are  to  be  subsequently  constructed  as  given  in 
Table  (8),  using  the  T's  and  the  two  blocks  of  D so  defined.  The 
appropriate  n^^-order  virtual-occupied  block  of  D is  obtained  by 
solving  Equation  (33)  using  the  resultant  X intermediate.  Finally, 
the  appropriate  nth-order  AO  derivative  part  is  obtained  by  using  the 
same  set  of  T intermediates. 

The  expressions  in  the  MBPT  tables  (Tables  (9-12))  are  in  terms  of 

ab 

the  low-order  iterates  of  T and  the  function  Q^(A,  BxC).  We 
summarize  the  low-order  iterates  of  A and  T and  define  related 
quantities  in  Table  (13). 

Summary  of  Computational  Strategy 

The  overall  strategy  for  computing  MBPT/CC  gradients  is: 

1)  Obtain  SX  and  QX  for  all  3N  degrees  of  freedom.  (See  Equations 

(25)  and  (26).) 

2a)  If  an  nth-order  MBPT  gradient  is  desired,  solve  the  unperturbed 

MBPT  problem  by  low-order  iteration  of  the  CC  t-amplitude 
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Table  9.  Summary  for  the  MBPT(2)  Model. 


“ij  ' 2 'M(1> 

Dab  * - 2 «“<»  'iS'1* 

r(ij , ab)  . j t’J(l)  i 

Note:  Summation  over  all  but 


all  other  T's  = 0 
target  labels  is  implied. 


Table  10.  Summary  for  the  MBPT(3)  Model. 


D.  . = 

ij 

2 'ki<2>  * 

1 

2 

“jk<2>  'M*1* 

D . = - 

ab 

2 ‘S'”  *i3b<2>  - 

1 

2 

At*J<2)  tjJ(l) 

T(ij ,ab) 

- 2 ^ 

T(ab, cd) 

- B tij(1) 

r(kl,ij) 

= B 

*S(‘> 

T(ja,bi) 

t ik.( 1 ) 

5 

T(ai , be)  = I\jk,ia) 

Note:  Summation  over  all 

but  target  labels  is  implied. 

= 0 
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Table  11.  Summary  for  the  MBPT(4)  Model. 


DU  - - 

t j (2)  tf(2) 

1 

+ 2 

t^d)  tg(3)  + 

1 

" 12 

< jlk<23  ‘lki<2) 

D . = 

ab 

t*(2)  tj(2) 

1 

" 2 

t*J(l)  tjJ(3)  - 

1 

+ 12 

tijk(2)  ^k<2) 

2 «?k<2>  'kj<2>  + 2 64b(3)  ‘ki™ 


1 stac 


cb . 


2 -i3<2>  - 5 CXac<3>  'u(1) 


.cb. 


r(ij  »ab) 

= 

1 

2 

<u‘3> 

1 

+ 5 

Qjj(t(l),  t(l)xl)  / 

1 

+ 5 

ab 

(ten, 

ixt(i))/  : 

Dij  + 5 Qij(I’  t<1>xt(1)) 

* l 

Qjj(t(l), 

t(l)xt(l)) 

T(ab, cd) 

= 

1 

8 

ta5<i) 

tfj (2)  + 

8 StiJ(2)  'ij(1) 

r(kl,ij) 

= 

1 

8 

t^o) 

tg(2)  ♦ 

8 4t?j<2>  'n*13 

T(ja,bi) 

= 

‘w<2> 

. A«“<2>  tg(l) 

T(ai,bc) 

= 

ta(2) 

*s<*> 

* 2 'ik3<2>  'k?<13 

T(jk,ia) 

= 

- 

tj(2) 

*>> 

- 2 '3k5‘2>  ‘h'1' 

Note:  Summation  over  all  but  target  labels  is  implied. 
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Table  12.  Summary  for  the  LCCD  (or  D-MBPT(°°))  Model. 


D.  . = 

ij 

1 

2 

*£< 

“) 

D , = - 

ab 

1 

2 

*«< 

•>  ‘3° 

") 

T(ij ,ab) 

= 

1 

2 

*«<•> 

T(ab, cd) 

= 

1 

8 

*3 

r(kl,ij) 

= 

1 

8 

*3<-> 

*S 

T(ja,bi) 

= 

*£<-> 

*3 

r(ai,bc)  = r(jk,ia)  = 0 


Note:  Summation  o 

The  symbol  t?j 


ver  all  but  target  labels  is  implied. 

(oo)  represents  the  converged  LCCD  amplitudes. 
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Table  13.  The  Low-order  Iterates  of  A and  T. 


tij(i>  = xib(i>  = <ijMab>/Di$ 

D^b  = e.+  e.-e-e. 

ij  i j a b 


t??(2)  = X^(2)  = tJJ(l)  + At®?(2) 


ab. 


tab. 


'ij 


ab 


iJ 


ij 


AtjJ(2>  = [l  tjj(l)  <ab | | cd>  + \ tJJ(l)  <i j | | kl> 

+ P(ab|ij)  tg(l)  <jc|  |bk>  ] / DjJ 


tij(3)  = tij(1)  + Atij(2)  + AtiJ(35SDT)  + AtjJ(3;Q) 

xab(3)  = tfj(l)  + AtjJ(2)  + At^(3;SDT)  + AX*j(3;Q) 


Lij 


ij 


AtjJ(3;Q)  = QjJ(I,  t(l)xt(l))  / DjJ 


AX*j*(3;Q)  = {Q^(t(l),  Ixt(l))  + QjJ(t<l),  t(l)xl)}  / ojJ 


t(l)  = tjJ(l)  1=  <ij I I ab> 


-,ab. 


^bd  -bd  „ac 


cu  - 2 1 Blj  ci:i  * »« cki ) 


Q““(A,  BxC)  - \ A'J  [ eg  - 2 ( B« 

- 2 ( Bik  Cjl  + *fl  Cik  ) * 4 ( Kl  C“  * Bki  Clj 

At^(3;SDT)  = [j  At?j(2)  <ab||cd>  + j AtjJ(2)  <:ijllkl> 


+ P(ab | ij  ) At“(2)  <jc|  |bk> 


+ P(ij)  t?(2)  <cj | | ab>  - P(ab)  tj(2)  <ij||kb> 

\ P(ab)  tj^(2)  <cd  | | bk>  - \ P(ab)  tgj(2)  <jc|  |kl>  ] / V*h. 


Table  13 — continued 


t*(2)  = 

X*(2)  - 

[-!<>  <b£' 

|aj>  + \ tjj(l)  <ib|  |jk> 

D?  - 

z.-  e 

1 

i a 

^bk<2> 

■ Xabc<2> 

= [ - P(a/bc|k/ij)  tj*(l)  <dk| |bc> 

♦ P(c/ab| i/jk)  tjJ(l)  <jk| | lc>  ] / 


Dijk  " V V V V V Ec  • 


Note:  Summation  over  all  but  target  labels  is  implied. 
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t h 

equations.  Save  the  t-amplitude  iterates  through  n-1  order  in 
the  correlation  perturbation. 

b)  Or  if  a CC  gradient  is  desired,  solve  the  unperturbed  CC  problem 
by  converging  the  CC  t-amplitude  equations.  Save  the  converged  t- 
amplitudes . 

3a)  If  an  MBPT  gradient  is  desired,  construct  the  T intermediates  and 
the  occupied-occupied  and  virtual-virtual  blocks  of  the  density 
matrix  D from  the  low-order  iterates  of  the  CC  t-amplitude 
equations.  (See  Tables  (9-12).) 

3b)  Or  if  a CC  gradient  is  desired,  it  is  necessary  to  solve  the 
linear  equation  for  the  X-amplitudes  using  the  converged  t- 
amplitudes.  Then  construct  the  T intermediates  as  given  in  Table 
(6)  and  construct  the  occupied-occupied  and  virtual-virtual  blocks 
of  the  density  matrix  D as  given  in  Table  (5). 

4)  Transform  the  T intermediates  to  the  AO  basis;  compute  the  AO 
derivative  contributions  for  each  of  the  3N  degrees  of  freedom. 

(See  Table  (7) . ) 

5)  Using  the  T intermediates  and  the  occupied-occupied  and  virtual- 
virtual  blocks  of  the  relaxed  density  matrix  D,  construct  the 
intermediates  I and  X as  given  in  Table  (8). 

6)  Solve  Equation  (33)  to  obtain  the  virtual-occupied  block  of  the 
density  matrix. 

7)  Evaluate  the  simple  dot-products  in  Equation  (34)  for  each  of  the 
3N  degrees  of  freedom;  add  to  the  corresponding  AO  derivative 
contributions  to  obtain  the  total  derivative  of  the  correlation 
energy  for  each  degree  of  freedom. 

v 

(Notes  In  Equation  (34)  the  density  D,  the  perturbation  matrix  Q 
and  overlap  derivative  matrix  S^  are  in  terms  of  the  MO  basis.  The 
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3N  sets  of  SX  and  QX  matrices  are  initially  prepared  in  the  AO 
basis  before  an  MO  transformation.  (See  Equations  (25)  and  (26).) 
It  is  more  efficient  to  avoid  3N  sets  of  transformations  by 
instead  transforming  the  perturbation-independent  matrices  D and  I 
to  the  AO  basis.  The  contributions  of  Equation  (34)  can  then  be 
computed  equivalently  as  dot  products  of  matrices  in  the  AO 
basis. ) 

8)  Option:  Since  the  MBPT/CC  relaxed  density  is  available,  one  may 
readily  compute  one-electron  response  properties  or  plot  the 
electronic  density  of  the  molecular  system. 


SECOND  DERIVATIVES  OF  THE  COUPLED-CLUSTER  ENERGY 


Derivation  of  General  Expressions 


Let  us  again  return  to  Equation  (6),  where  the  Hamiltonian  is  in 
normal-ordered  form  and  the  wavefunction  is  the  exponential  of  the 
Coupled-Cluster  model.  Let  there  be  two  perturbation  parameters  a and 
P,  upon  which  the  Hamiltonian,  the  wavefunction,  and  AE  are  dependent. 


;-T(a,  0)  eT(a’P)|0>  = AE(a,P)|0> 


(35) 


Let  us  define  the  following: 


T = T(a,  P) 


„ 3T(a,  P) 

Ta 

3T(a,  P) 
TP 

a=P=0 

T 

3a 

a=P=0 

90 

e 32K“.S) 

1 - ' 

3a3p 

a=  P=0 

a=  0=0 


HN  = fN  + WN  = fN^a’^ 


a=p=0 


+ VN(a»P) 


a=P=0 


Bj  - £“  + V“  . 


f9fN(a,p)  8VN(a,p)) 
+ 


k 3a  3a 


a=p=0 


hn  - tN  N " 


f3fN(a,p)  3VN(a,p)^ 
+ 


3P  30 


a=  P=0 
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ap  _ fap  „aP 
HN  " £n  wn 


p2fN(a,P)  32WN(a,P)-| 


, 3a3P 


3a3p 


a=P=0 


AE  = AE(a, P) 


AEa  = 


3AE( a,  P) 


06=  3=0 


3a 


AEP  = 


a=  P=0 


3AE(a,P) 

3P 


a=p=0 


AE 


ap 


3^AE(a,p) 

3a3p 


«x=P=0 


As  before,  projection  on  the  left  of  Equation  (35)  by  <P|  and 
evaluation  at  <x=P=0,  yields  the  unperturbed  Coupled-Cluster 
expressions.  The  derivative  of  Equation  (35)  with  respect  to  a (or 
P),  evaluated  at  ot=P=0,  is  the  equivalent  of  Equation  (14)  of  the 
previous  section: 


KhN  eT>c’  + (HN  *T)c 


|0> 


Or  more  simply, 


AEa|0>  . 


(36) 


|0> 


AEa|0>  . 


Projection  on  the  left  by  <P|  yields  the  equivalents  of  Equations  (15) 
and  (16):  an  expression  for  AEa  and  a set  of  linear  equations  which 
determine  the  derivative  t-amplitudes: 


<0|((Hn  eT)cTa)c|0>  + <0|(H“  eT)c|0>  = AEa  (37) 


and 

<*!((%  eT)cTa)c|0>  + <*|(Hj  eT)c|0>  = 0 . (38) 
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Although  the  evaluation  of  first  derivatives  does  not  require  the 
derivative  t-amplitudes,  the  evaluation  of  second-derivatives  does. 
The  coupled,  linear  equations  represented  by  Equation  (38)  for  the 
CCSDT-1  model  are  presented  in  Tables  (14-16). 

The  second  derivative  of  Equation  (35)  with  respect  to  a and  0, 
evaluated  at  a=S=0,  is  given  by: 


[(HNeT)c,  Taf5]  + l(H^eT)c,  Ta]  + [(H“eT)c,  T*5] 


(39) 


+ (HNPeT)c  + tt<HNeT)c’  TtX]’  TP) 


1 0>  = 


AEag|0> 


Thus , 


«H„0.  T“S>„  ♦ ((H“e‘)„  T“)„  * T") 


ST,  -a. 


,a  T , 


N 'c 


N ' c 


N ' c 


Tx 

(HN  e >c  + 


(((HNeT)cTa)cTg)c 


1 0>  = 


AEaP|0> 


Projection  on  the  left  by  <P|  will  yield:  an  expression  for  AE  p in 
terms  of  Ta,  T^,  the  unperturbed  T,  and  Ta^;  and  a set  of  coupled 
linear  equations  for  the  Ta^  amplitudes.  The  dependence  of  AEa^  on 
Ta^  can  be  eliminated  just  as  the  first  derivative  T is  eliminated 
from  the  first  derivative  of  the  energy.  That  is,  we  return  to 
Equation  (39)  and  insert  1 = |PXP|  into  the  terms  containing  T . 
Since  <0|TaP|P>=0  and  < | ( HNeT ) c 1 0>=0 , we  obtain  upon  left-projection 
by  | P> : 
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<0| 


«9{jeT)cTa)c  * ((Hje‘)c  T°)c  * <H*>V)c  * 


,a  Tv  t3v  . ,na6  T 


Tv  rr<Xv  TPv 


|0> 


+ <0|(HNeT)c|*X$|TaP|0>  = AEag 


and 


<$|-TaP|0>  AE  + <*|(HNe1)c|$X$|T“M|0> 


ag. 


+ <*| 


( ( HNgT > c T“>c  + ((HNeT)c  T*>c 


[A] 


IB] 


+ wf*T)c  + (((HNeT)cTa)cTe)c 


|0>  = 0 


[C] 


ID] 


a8 

The  second  equation  is  solved  for  <f|T  p|0>: 

<*|TaP|0>  = - <*|  (HNeT)c-  AE  | *>  <*| ( [A]+[B]+[C]+[D] ) |0>  , 

and  we  substitute  for  <*|Ta^|0>  in  the  AEa^  expression: 


<0| 


(<H^T)ct<I)c  - (("^T>c  T\  + <HNSeT)c  + <«V‘>CT">CTH) 


Tv  fpOtv  T0v 


|0> 


-1 


- <0 | ( (HNeA )c | *>  <*| (Hn  ei)c~  AE  | *>  <i | ( [A]+[B]+[C]+[D] ) | 0>  = AE 


a3 


Recalling  the  definition  of  the  perturbation-independent  de-excitation 
operator  A,  we  have 
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\ 


\ 


+ <0|A|*>  <*|([A]+[B]+[C]+[D])|0> 


AE 


aP 


An  equivalent  form  of  the  above  expression  can  be  readily 
subjected  to  the  usual  diagrammatic  rules.  Recalling  that  1 = |PXP| 

= |0><0|  + |*X*|,  the  second  derivative  of  the  energy  can  be  written 

as 

AEaP  = <0|((H“eT)cTP)c|0>  + <0( (HPeT)cTa)c | 0>  (40) 

+ <0|(H“PeT)c|0>  + <0|(((HNeT)cTa)cTP)c|0> 

+ <0| A ((H“eT)cTP)c|0>  + <0|A  ((HjeT)cTa)c|0> 

+ <0|A  (HjfeT)c|0>  + <0|A  (((HNeT)cTa)cTP)c|0>  . 


Equation  (40)  is  the  general  expression  for  the  second  derivative 
of  the  CC  energy.  It  applies  to  CC  models  of  any  truncation  of  T or 
choice  of  spin-orbital  basis.  We  wish  to  develop  the  general  second 
derivative  energy  expression  for  the  specfic  choice  of  the  CCSDT-1 
model.  The  first,  second,  and  fourth  terms  of  Equation  (40)  are 
easily  evaluated,  as  they  reduce  to  forms  analogous  to  the  ordinary  CC 
energy  expression: 
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1 

5 

<ab  | 

1 ij>“ 

ab 

+ 

<ab  | 

|ij>“ 

1 * 

f.  axa 
ia  i 

1 

+ 4 

<ab  | 

|ij>e 

Tab  + 

1J 

<ab  | 

|ij>P 

*: 

J 

f.  6 T» 
ia  i 

+ 

<ab | | ij> 

< 

xb 

J 

where  the  symbol  x represents  the  derivative  of  a t-amplitude  with 
respect  to  a,  evaluated  at  a=8=0  and  x represents  the  derivative  of  a 
t-amplitude  with  respect  to  8>  evaluated  at  a=3=0. 

The  third  and  seventh  terms  of  Equation  (40)  are  exactly  the 
expressions  in  Table  (1),  except  that  the  symbol  X is  to  be 
interpreted  as  the  derivative  with  respect  to  a and  3,  evaluated  at 
cu=3=0.  Of  course,  Table  (1)  is  for  the  CCSDT  model,  but  we  may  simply 
ignore  terms  in  Table  (1)  that  do  not  apply  to  the  truncated  CCSDT-1 
model.  The  algebraic  expressions  for  the  sixth  term  in  the  CCSDT-1 
model  are  given  in  Table  (17).  The  fifth  term  is  obtained  from  the 
sixth  term  by  exchanging  a and  3*  Finally,  the  eighth  term  is  given 
in  Table  (18) . 

Expressions  for  the  Choice  of  SCF  Reference 

As  in  the  section  on  first  derivatives  of  SCF  orbitals,  all 
orbitals  are  assumed  to  be  unperturbed,  i.e.  evaluated  at  a=3=0, 
unless  otherwise  shown  as  a function  of  a and  8-  Derivative  orbitals, 
evaluated  at  X=0,  are  indicated  by  the  superscript  a8- 


1 P(a,  8)>  - 4»p(a,8) 

|p>  = |p(a,e)>la=p=0 

Ip“p>  - 

82|p(oc,  8)> 

3a38 

ot=  3=0 

|y(a,  8)>  = ^y(a»3) 

|y>  = |w(a,P)>la=p=0 

|ga6>  - 

a2lu(a,8)> 

3a38 

a=3=0 

The  results  of  CPHF  theory  for  the  second  derivative  of  non-canonical 
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Table  14.  The  T^a  Equation  for  the  CCSDT-1  Model. 


0 = 


f.  a + £.  a tb  - f.. 

ia  ba  l ji 

f a ..ba  .b  .a. 

+ fjb  (tji  - li  tj) 


“ tj  + <aj | | ib>a  tb 


\ <ja|  |bc>“  Cb?  - * <jk| |bi>a  Cb*  - 2 <jk|  |bc>“  t 
2 <jk||bc>a  tj  tjj  - <jk||cb>a  t£  (tfj  - t?  tb) 

l <bc||Jk>«  tj“ 


l,  i.  - j...  T®  + <aj||ib>  t|? 

ba  l j i j J 11  J 


S - £j 

ba  _b 


*;-*5 

$ * 

tc  + tb 

T? 

1 

- T? 

J 

t5  - t° 

’J) 

ti  + 1. 

ti  t. 

.ba  Tb 
jk  J 

..a  4.b 

*k  * ‘j 

i 

*:  - *j 

<> 

.b  tac 
i \Jk 

t*  t“> 
i jk 

.a  tbc 

■k  'ji  + 

<*s  - 

*j  *i>  * 

*£ 

(tab 

1J 

- t 

a b, 


+ 5 <bc | | jk> 

where  we  define 


= tab  ♦ t?  tb  - tb  t? 
1J  1J  1 J 1 1 


b ac 
i jk 


Note: 


Summation  over  all  but  target  labels  is  implied. 
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Table  15.  The  T2“  Equation  for  the  CCSDT-1  Model. 

0 . <ij||at»“  t (£bdet  + <bk|  |dc>“  t£) 

t (fM“+  <ak||dc>“  tj)  t?* 

- <fbd“.  <lk||cj>“tj)  «£ 

- <£bd“+  <lk||ci>“  'l>  'kj 


t? 

k 1J 


>“  cjj 


ad 

d 

. ab  td 

ab 

‘ij 

+ *i 

lkj  + 

lik 

+ 

<ij||kl>“ 

,ab 

kl 

* ‘5 

<kl | |ic>“ 

rab 

0kl 

- <bk|  |jd>a  cjJ  + <bk| |id>“  cjJ  + <ak|  |jd>“  - <ak|  |id>“  C^b 


"kj  T 'CW'-|IJV"  -kj 
lad  .c  /L,.i  i _ Jv  a ^,ad 


'jk 


- t?  <bk| |cd>“  C^j  * t'  <bk||cd>“  c“ 

- t?  <ak|  [cd>ot  + t?  <ak|]cd>“ 

. t»  <lk||cj>“  tjjj  - tj  <lk||ci>“  t‘  <bl||cd>*E^ 

- tj  <lk||cj>“  t“  . tj  <lk||ci>“  t“  -l  t*  <al||cd>*E^ 

- tb  <ak| |ij>“  + t®<kbf|ij>a+  t?<ab||cj>“+  t?<ab||ic>a 

1 y\r"\  I J r p^^  p^b  j , 3.C  pbd  *bd  Cac) 

+ 5 <kl  I I cd>  [ C.J  Ckl  - 2(ti;j  Ckl  + t..  Ckl) 

-*(‘S  <fi  * ‘ji  cik>  * 4<cK  <$$♦  <i  CU> 

-8  tf  t*  tj  1 

2 P(a/b)  <cd  | | bk>a  t£j“  - j P(i/j>  <1c||kl>“  t“j 


, a „cab 
* £ck  'kij  * 


r ad  ,,,  , i , v ,.c  aa  .c  au. 

+ fbd  Tij  + <bk||dc>  (tk  Tij  + lk  Tij} 

, db  .lit.  v ,.c  db  .c  dbv 

♦ £bd  Tij  * <ak| |dc>  <tk  Ty  * tk  Tdj) 

_ ab  i i ..  /.c  ab  .c  ab. 

‘ fbd  Tik  " <lk I I CJ>  <li  Tik  + Tik) 

r sb  yT  i it  »v  / . c ab  ,c  ab  v 

“ fbd  Tkj  + <lkHcl>  (tl  Tkj  + tl  Tkj) 


c ad 


.c  _ad. 
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Table  15 — continued 


<< 

*5* 

4 

tad  + 

ij 

Td 

1 

ab 

lkj  + 

Td 

J ik 

«k 

Tdb* 

1J 

4 

Tad  + 

IJ 

td 

1 

+ 

d ab 
t7  TV, 
J ik 

kd 

Ekd 


<ab|  |cd>  (C?d)a  + <ij||kl>  (CaJ)a 


xj  <kl||cj>  + xj  <kl||ic>  eg 


+ tj  <kl|  |cj>  (C^)“  + tj  <kl | | ic>  (Cg)a 
- <bk||jd>  (Cad)a+  <bk| | id>  (Cad)a  + <ak| | jd>  (Cd£)“  - <ak||id>  (cJJ) 


<bk  | | cd>  Cjj  «■  T?  <bk  | | cd>  Cad 
cj<ak||cd>  Cd£ 


x?<ak||cd>  Cd|j  + 


t?  <bk|  |cd>  (Cad)“  + tj  <bk||cd>  (C^)a 
t?  <ak|  | cd>  (Cjj)“  + tj  <ak||cd>  (Cd£)“ 


<lk| | c j > t?k 
<lk| |cj>  tJJ 

tj  <lk| | cj > 


x*  <lk|  |ci>  + J t*  <bl  | | cd> 


^ <lk|  | ci>  - l 
cb 


<al | | cd>  Ecd 


ij 

,cd>  a 


ta  <lk||ci>  t"  .J  tj  <bl | | cd>  (E?°) 

- tJ<lk||cj>T«  * tb  <lk||ci>  t“  -\  tb  <al|  |cd>  (E^d 

- x|j  <ak|  |ij>  + Ta<kb||ij>  + x?  <ab  | | c j > + T^<ab||ic> 

1 ..  . i i r ,_cd,.a  „ab  0,  ac  rbd  bd  rac, 

+ 5 <kl | | cd>  [ (C..)  Ckl  - 2(t..  Ckl  + T..  Cfcl) 


cd  ,cab.a 
+ C.j  (oklJ 


2(,-  (c“)“  ♦ «£)»> 


ab  rcd  — ^b  pCd v pbd  ._bdv<x  p3c.. 

’2(Tik  Cjl  + Tjl  Cik)  + 4((Cki)  Clj+  (Cki>  Clj} 


-2(t^  «$>“  ♦ tab  <>“)  * MC-  (Cbd>«  * Cbd  «")“) 

M «k  *i  ‘J  * ‘J  -k  *i  '5  * 'i  4 4 >5  * 4 4 «*  *5» 

£ck  Tkij  * j P<a/b>  <cd||bk>  Tkij  - 2 P(i/J)  <Jc||kl>  Tkli 


103 


Table  15 — continued 


where  we  define 


rab 

I35  + 

ta 

tb 

ta 

cij  ■ 

tij  + 

li 

J 

= Tab 

1J 

x 

Ta 

1 

tb  _ 

Tb 

1 

T 

J 

lab 
Cij  = 

<5 

_ab 

Eij 

*S- 

‘1 

t® 

3 

(E“b)“ 

1J 

- <5 

- 

tS + 

Tb 

1 

ta 

3 


ta 

l 


rb 

3 


l 


ca 

3 


(ca5>«  = 


i 


$ 


ab  a b 
Tij  + Ti  *3 


'i’J 


Note: 


Summation  over  all  but  target  labels  is  implied. 
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Table  16.  The  T3“  Equation  for  the  CCSDT-1  Model. 

0 - - P(k/iJ)  tjjj  fkl“  - P(k/ij)  tjy  fkl 
+ P(c/ab)  £dc“  ♦ P(c/ab)  £dc 

- P(a/bc|k/ij)  tf?  <dk||bc>“  - P(a/bc|k/ij)  Tj?  <dk||bc> 

1 j J 

+ P(c/ab | i/jk)  tjJ  <jk||lc>a  + P(c/ab|i/jk)  <jk||lc> 


Note:  Summation  over  all  but  target  labels  is  implied. 
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Table  17, 


+ 


+ 


+ 


+ 


+ 


+ 


Contributions  to  the  Second  Derivative  of  the  CCSDT-1 
Energy  from:  „ T 

<0 | A ((Hjei)cTa)c|0>. 


Xh  fh  ^ 
b ba 

V*  f..6 

a ij 

<ib|  |aj>^ 

xj  «b 


3 a 


Xj  f e tb 
a £ib 


X^  <ci 
c 

x*  <CJ 

Xa  ^ 
Xb  <ki 

K <Jk 

Xj  <ki 


XJ  <ki 
c 

xc  <k3 

Xk  <ji 
a J 

<cd 

cd 

xkl  ... 
Xba  <J1 

xki  . . 
Xcb  Oc 

Xkj  <ic 

ca 

Xbl  <bc 


ba>e  tb 
ba>*  tb 


ka>0  tb 


kb>P  £j 
N.P 

Ca>  'jk 

bo^  $ 

ba>^  t£  t 
ba>^  t£  t 

cb>B  *k  ' 
|ba>»  «J 

|kl>P  tb 

|ak>P  tj 

|bk>P  tb 
iaj>3 


b 

j 

b 

j 

b 

j 
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Table  17 — continued 


+ 

1 

2 

xjk 

ba 

<ib | |jk>P 

- 

1 

2 

xj1 

be 

fka* 

Jk 

1 

f P 

be 

2 

Xba 

ic 

'jk 

- 

1 

2 

Xjk 

Adb 

<di 

ca>^ 

be 

'jk 

+ 

Xji 

Adb 

<dk 

ca>^ 

tb^ 

jk 

- 

1 

5 

xjk 

da 

<di 

bc>P 

be 

jk 

+ 

1 

2 

be 

<ki 

la>* 

be 

'jk 

+ 

1 

5 

£ 

<jk 

la>P 

be 

'jk 

- 

<ki 

lO* 

be 

'jk 

- 

<dk 

ba>P 

c b 

*k  '3 

- 

1 

2 

<di 

cb>P 

*£ 

+ 

1 

2 

x1* 

cb 

<kj 

la>P 

*£  *S 

+ 

X« 

<ki 

lb>P 

t?  tk 

ca 

k J 

- 

1 

2 

Xkj 

cd 

<li 

ba>P 

cd  b 
'kl  J 

+ 

1 

2 

XkH 

cd 

<lj 

ba>P 

cd  b 
'kl 

+ 

1 

5 

xj* 

cd 

<kl 

ba>0 

cd  b 
'kl  J 

- 

1 

2 

xki 

cb 

<ji 

da>P 

cd  b 
'kl  'j 

+ 

1 

2 

xkl 

ca 

<ji 

db>P 

cd  b 
'kl 

- 

kcb 

<lj 

da>P 

cd  b 
'kl  lj 

- 

Xkj 

ca 

<li 

db>P 

cd  b 
'kl 

+ 

1 

5 

Xba 

<ji 

cd>P 

cd  b 
'kl  rj 
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Table  17 — continued 


1 

+ 2 

£ 

<lk| |ba>P  t*  t£  tj 

1 

+ 2 

<li| |cb>P  tj  tj  tj 

ab  [ 

+ T.  . X 

1 3 L 

Xs  £..6 
a jb 

+ 

X*  <kj  | | cb>P  tj 

1 

+ 2 

X*  <cj ||ab>P 

1 

" 2 

Xk  <ij||kb>P 

1 

" 2 

X*  <kj  | | ab>P  t 

1 

‘ 2 

Xk  <ij | | cb>P  t 

1 

+ 2 

Xk‘  f..6 
ab  ]k 

1 

“ 2 

Xij  f p 
ca  cb 

1 

+ 8 

X^  <cd  | | ab>P 

1 

+ 8 

>S 

+ 

Xk*  <jc| |bk>P 

1 

+ 2 

X«  £uk6  t? 
ca  kb  k 

1 

+ 2 

Xk*  £.  6 t? 
ab  jc  k 

1 

' 5 

Xjj  <dk| |ab>P 

+ 

X$*  <dj | |cb>P 

1 

+ 2 

Xjj  <dk| |cb>P 

- 

X^  <kj||lb>P 

1 

+ 5 

Xkk  <ij||lc>6 

1 

' 2 

X“  <kj||lc>B 

?r  n n 7?  <"> 
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Table  17 — continued 


- 

1 

4 

Xk* 

cd 

<lj| 

| ab>3 

cd 

rkl 

+ 

1 

15 

X1^ 

cd 

<kl  | 

| ab>3 

cd 

Tkl 

- 

1 

4 

Xkl 

ca 

<ij| 

|db>3 

cd 

£kl 

+ 

Xki 

ca 

<lj| 

|db>3 

cd 

Tkl 

- 

1 

4 

x1* 

ca 

<kl  | 

| db>3 

cd 

rkl 

+ 

1 

15 

ab 

<ij| 

| cd>  3 

cd 

£kl 

- 

1 

4 

Xkf 

ab 

<lj| 

| cd>3 

4t 

+ 

1 

8 

x*j 

dc 

<lk  | 

| ab>3 

4 

4 

- 

<lj| 

| cb>3 

4 

4 

- 

1 

2 

<lk| 

| cb>3 

4 

4 

+ 

1 

8 

<ij| 

| dc>3 

4 

4 

- 

1 

2 

<kj| 

| dc>3 

4 

4 

+ 

1 

4 

»51 <cd  i [bk> 

3 

- 

1 

4 

,kli  .. 
Xcab  <>' 

1 1 kl> 

6 ] 

abc 

ijk 

X 

[ \ 

Xk  <ij 1 |ab>3 

+ 

1 

4 

41 

f . 3 
kc 

- 

1 

4 

<dk| | bc> 

3 

+ 

1 

4 

xl£ 

ab 

<jk| | lc> 

3 

1 

xlij 

f 3 

12 

abc 

rkl 

1 

\ijk 

f 3 1 

12 

aab 

ac  J 

Note:  Summation  over  all  labels  is  implied. 

The  symbol  T represents  derivative  t-amplitudes;  the  derivative  is 
with  respect  to  a,  evaluated  at  a=0=O. 
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Table  18. 


Contributions  to  the  Second  Derivative  of  the  CCSDT-1 


Energy  from: 


<0 1 A (((HNeT)cTa)cTP)c|0>. 


l 


* [ 

- 4fja 

<5 

- 

xi  £ib 

Tb 

J 

+ 

<ci 

c 

| ba> 

<5 

- 

<cj 

| ba> 

- 

% <3i 

| ka> 

+ 

Xa  ^ 

|kb> 

Tb 

1 

+ 

Xj  <ki 

| ca> 

T*? 

jk 

1 

+ 2 

| ca> 

be 

Tjk 

1 

+ 2 

Xj  <ki 
a 

| bc> 

Tbc 

Tjk 

- 

Xj  <ki 
c 

| ba> 

'^b 

+ 

Xc  <kj 

| ba> 

+ 

.k  ... 
Xa  <J1 

| cb> 

*£^ 

- 

Xj  <ki 
c 

| ba> 

'X 

+ 

Xc  <kj 

| ba> 

'i* 

+ 

Xk  <ji 
a J 

| cb> 

tb  TC 
k 

1 

+ 2 

Xjj  <cd | | ba> 

Tb 

J 

1 

+ 2 

<^11^ 

Tb 

J 

X^<jc||ak>  Tb 
Xca<ic"bk>  Tj 


xj1 

be 


ka 


1 

2 
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Table  18 — continued 


- 

1 

2 

£ 

fic 

tbc 

jk 

- 

1 

2 

$ 

<di 

ca> 

T ^ 
jk 

+ 

xJi 

ab 

<dk 

ca> 

be 

Tjk 

- 

1 

3 

xjk 

da 

<di 

bc> 

xb^ 

jk 

+ 

1 

2 

£ 

<ki 

la> 

be 

Tjk 

+ 

1 

3 

<jk 

la> 

be 

jk 

- 

£ 

<ki 

lc> 

£ 

- 

xj1 

dc 

<dk 

ba> 

- 

xkj 

a a 

<di 

cb> 

+ 

X1' 

cb 

<kj 

la> 

c b 
t,  T . 

k j 

+ 

X1* 

ca 

<ki 

lb> 

- 

xj1 

dc 

<dk 

ba> 

ti 

+ 

X1* 

ca 

<ki 

lb> 

- 

1 

2 

Xkj 

cd 

<li 

ba> 

tcd  Tb 
*kl  Tj 

+ 

1 

2 

XkH 

cd 

<lj 

ba> 

cd  b 
*kl  Tj 

+ 

1 

3 

Xjj 

cd 

<kl 

ba> 

cd  b 
kl  Tj 

- 

1 

2 

XkJ 

cb 

<ji 

da> 

tcd  tb 

lkl  Tj 

+ 

1 

2 

Xkl 

ca 

<ji 

db> 

tCd  Tb 

*kl  Tj 

- 

xki 

cb 

<lj 

da> 

tCd  Tb 

*kl  Tj 

- 

Xkj 

ca 

<li 

db> 

tCd  Tb 
lkl  Tj 

+ 

1 

3 

£ 

<ji 

cd> 

tCd  Tb 

'kl  Tj 
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Table  18 — continued 


- 

1 

2 

»S 

•H 

rH 

V 

| ba> 

cd  b 
Tkl  3 

+ 

1 

2 

£ 

<lj 

| ba> 

cd  b 
Tkl 

+ 

1 

5 

cd 

<kl 

| ba> 

cd  b 
Tkl 

- 

1 

2 

cb 

<ji 

| da> 

cd  b 
Tkl 

+ 

1 

2 

Xkl 

ca 

<ji 

| db> 

cd  b 
Tkl  3 

- 

xki 

cb 

<lj 

| da> 

cd  b 
Tkl  3 

- 

Xkj 

ca 

<li 

| db> 

cd  b 
Tkl  3 

+ 

1 

5 

Xkl 

ba 

<ji 

| cd> 

cd  b 
Tkl  J 

+ 

1 

2 

xj1 

dc 

<lk 

| ba> 

td  tc  Tb 

rl  rk  Tj 

+ 

Xkj 

da 

<li 

| cb> 

td  tc  tb 
rl  rk  Tj 

+ 

xj1 

dc 

<lk 

| ba> 

Td  tC  tb 

T1  rk  rj 

+ 

1 

2 

XHj 

da 

<li 

| cb> 

xd  tc  tb 
T1  lk 

ab 
+ -“f.  . 

ij 

* I 

X*  <kj | | cb> 

- 

1 

2 

X*  <kj | | ab>  T 

- 

1 

2 

Xk  <ij | |cb>  t 

+ 

1 

2 

Xij 

ca 

f.  . X? 
kb  k 

+ 

1 

2 

xki 

ab 

be  Tk 

- 

1 

5 

X*j 

dc 

<dk| | ab> 

+ 

Xki 

da 

<dj | | cb> 

+ 

1 

2 

xjj  <dk||cb> 

X11 

ca 

<kj | | lb> 
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Table  18 — continued 


1 

+ 3 
1 

" 2 
1 

' 3 
1 

+ 15 
1 

" 3 


1 

‘ 3 
1 

+ 16 
1 

' 3 
1 

+ 3 


xlk  . , 

Xab  <1J  I 

X1,1  <kj  | 


1 

2 

1 

2 

1 

3 

1 

2 

1 

2 


ab 

xc5  i 

Xcd  <kl' 

Xkl  <i  J I 
ca  J ' 

Xki  <lj| 

Xij  <kl | 
ca 

Xa£  <»  I 

X"b  «JI 

Xdc  <lk' 

Xdl  <«' 
x£  <«  i 
xdi  <iri 
xJa  <lk| 

Xab  <1J I 
Xab  <kj I 
Xab  <k5  1 


|lc> 

|lc> 

| ab> 

| ab> 

| db> 
|db> 

| db> 

| cd> 

| cd> 

| ab> 

| cb> 

| cb> 

| cb> 

| cb> 
| dc> 
| dc> 
| dc> 


lk 

cd 

Tkl 

cd 

Tkl 

cd 

Tkl 

cd 

Tkl 

cd 

Tkl 

cd 

Tkl 

Tcd 

Tkl 

td  TC 

rl  Tk 

td  TC 

ll  Tk 

Td  tC 
T1  \ 

td  TC 

ll  \ 

Td  tC 

T1  zk 

td  TC 

rl  Tk 

td  TC 
Z1  Tk 


td  tc  1 

rl  Zk  J 


Note:  Summation  over  all  labels  is  implied. 

The  symbol  x represents  derivative  t-amplitudes;  the  derivative  is 
with  respect  to  a,  evaluated  at  ot=3=0. 

The  symbol  x represents  the  derivative  t-amplitudes  where  the 
derivative  is  with  respect  to  g,  evaluated  at  a=p=0. 
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SCF  orbitals  {<|>(a,  P)},  expanded  in  a primitive  AO  basis  set  {»Ha,  3)} 
are  now  summarized. 

We  follow  Handy  et  al.,  and  choose  non-canonical  orbitals  such 
that  the  second  derivative  of  an  orbital,  evaluated  at  oc=P=0,  is  given 
by  a linear  combination  of  unperturbed  SCF  orbitals  where  there  are  no 
singularities  arising  from  degeneracies  among  orbital  eigenvalues.  Of 
course,  there  are  contributions  from  the  second  derivative  of  the 
primitive  AO  functions,  as  well. 


|aoB> 

liae> 


l * l * J 

£ U“f|t>  ♦ E U“f|k>  + E c .|u“B> 

ft1  k y p 


where 


,“3 

Jfa 


1 nap  1 Fag 

- 2 ha  and  Uki  = "2  ^ki 


The  coefficients  UaP  are  the  second-order  CPHF  coefficients.  The 

pq 

virtual-occupied  block  of  the  CPHF  coefficients  are  determined  by 
solving  the  second-order  CPHF  equation: 


Uai  = 


aP 


Ej  " Eb 


where 


A . , . = 1 

ai,bj 


<ab I I ij>  + <aj 1 I ib> 


E - S. 
a l 


Note  that  the  matrix  A is  the  same  as  in  the  first-order  CPHF  equation 
given  in  Equation  (21).  The  occupied-virtual  second-order  CPHF 
coefficients  are  determined  by  the  second  derivative  of  the  SCF 
orthonormality  condition: 


UaP 

pq 


+ 


(uaV  + 

qp 


o , 
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where 

rag  S«P  £ y3  ytt  yB 

pq  pq  r Pr  Q r <lr  Pr 

fY.fi 

The  elements  of  the  matrix  Q p are  given  by: 


Sa  SP  - SP  Sa  ) 
pr  qr  pr  qr' 


cxg  hap  _ s«e  £ _ z S“P  <pl  | | qk> 

pq  pq  pq  pq  k,i  kl 


+ E u“  {<ps|  |qj>P  + <pj|  Iqs>3} 


sj 


S,j 

+ E uf.  { <ps  I |qj>a  + <pj  | |qs>a} 


s,  j 


sj 


+ I U“  U*.  { <ps | | q t > + <pt| |qs>} 

s,t,j  J J 

+ £ Ua  f P + EUa  f P + E Ua  f P + E U?n  f “ 
sp  sq  s sq  ps  g sp  sq  g sq  ps 


- E 
s 


[ Ua  Up  e + U*3  Ua  e + U*  uj3  e 

l sp  qs  q sp  qs  q sq  ps 


+ UH  IT  e + Ua  UP  e + U*  ec  1 

sp  sq  s sp  sq  s J 


jP  y« 

sq  ps  p 


- (c  p 3<y(«,3)X(«.3)l  Iv(«,  P)g(«,  3)> 

+ ' yp'  vq  Xcr  3a33 

y»v  * ^ 

X,  a 


a=(3=0 


where 


'pq 


1 . . when  p and  q are  both 

2 'ep  + virtual  or  occupied 

s otherwise 


The  matrix  ha^  consists  of  the  derivative  of  the  matrix  elements  of 

pq 

the  one-electron  part  of  the  Hamiltonian  in  the  primitive, 
perturbation-dependent  basis  {y(cc,3)},  but  transformed  into  the 
unperturbed  orbital  basis  {<(>}. 


haP  = E (C  )*C 

pq  y>v  yp  vq 


3<y  («,  P)  | h(a,  g)|v(a,  p)> 


3a33 


a=  3=0 


Similarly,  the  matrix  is  the  derivative  of  the  overlap  matrix  of 
the  perturbation-dependent  primitive  basis  functions,  transformed  into 
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the  unperturbed  orbital  basis. 


E 

-y,  v 


^Cyp^  C\>q 


9<u(g,  g)  | v(g,  g)> 


g=g=0 


For  the  choice  of  non-canonical  orbitals,  the  off-diagonal  elements 

f°f  and  f??  are  not  zero.  These  matrix  elements  are  given  by: 
ab  ij 


fa“P  = + 2 U°?(<am| |bg>  + <ag| |bm>) 

a g>m  8 

f“P  = Qj?  + 2 U°?(<im|  |jg>  + <ig||jm»  . 

J J g,m  g 

Recall  that  the  third  and  seventh  terms  of  Equation  (40)  are 
analogous  to  the  first  derivative  of  the  energy  presented  in  Table 
(1).  Now,  using  the  identities  from  CPHF  theory,  we  may  expand  the 
derivative  integrals  just  as  in  the  previous  section.  We  will  obtain, 
in  analogy  with  the  first  derivative,  the  following: 


A E 


gg 


l D.  . Q“?  * £ D . Q“f  ♦ 2 £ D , QT  (41) 

l,i  13  13  a,b  ab  ab  a,  i al  “ 


,ag 


ag 


+ E I.  . + 2 I , C + 2 E I,  ( 

i,3  J a,b  a,  l 


.gg 


gg 


+ A0/A0,  A0/M0  and  M0/M0  derivative  parts 

+ \ 2 <ab||ij>a^  + 2 <ab||ij>ata^ 

i » j » 1 » J * 

a,b  a,b 

+ l E <ab||ij>P  + 2 <ab||ij>P  ta  xb. 

i > j > 1 > j » 

a,b  a,b 

+ E <ab | | ij>  Ta  ^ 

i.j.  J 

a,  b 


+ Terms  5,6  and  8 of  Equation  (40) 
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The  intermediate  I and  the  density  matrix  D are  exactly  as  defined  in 

the  previous  section.  The  AO/AO,  AO/MO  and  MO/MO  derivative  parts, 

not  explicitly  shown  above,  are  presented  in  Table  (19).  Recall  that 

terms  5,  6 and  8 of  Equation  (40)  are  summarized  in  Tables  (17)  and 

(18)  and  that  the  derivative  t-amplitudes  are  determined  by  solving 

the  coupled,  linear  equations  given  in  Tables  (14-16).  Of  course,  the 

matrix  elements  f .,  f.  , f .a  and  f.  a are  zero  for  the  choice  of  SCF 
ai  la  al  la 

basis  functions,  wherever  they  may  occur  in  the  tables.  The  matrix 

elements  f ua  and  f . are  given  in  Equations  (30)  and  (31).  Finally, 
ab  ij 

the  derivative  integrals  which  appear  in  Equation  (41)  are  dependent 
upon  the  first-order  CPHF  coefficients.  The  derivative  of  a general 
integral  <pq||rs>  is  given  by: 


<pq I I rs> 


* * 3<y(a,  P)v(a,  P)  | p(a,  8)  cr(a,  g)> 


Z(CC  - C C ) C C 
v pr  as  ps  ar/  yp  vq 


3a 


y,  v, p, a 


a=P=0 


+ Z 
t 


a a a a 

<tq | | rs>  U + <pt | | rs>  U + <pq||ts>  U + <pq||rt>  U 

tp  tq  tr  ts  . 


MBPT  Second  Derivatives 


Now  we  extract  the  finite-order  MBPT  second  derivatives  from 
Equation  (41).  Inspection  of  the  first  six  terms  of  Equation  (41) 
reveals  that,  as  with  the  first  derivative,  the  second  derivative  of 
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the  nth-order  correlation  energy  depends  upon  the  nth-order  MBPT 
density  matrix  D and  the  nth-order  intermediate  I.  The  finite-order 
MBPT  densities  and  I intermediates  are  discussed  in  the  previous 
section;  Tables  (9-12)  prescribe  the  T intermediates  necessary  for 
their  construction.  The  appropriate  nth-order  AO/AO,  AO/MO  and  MO/MO 
derivative  parts  of  Equation  (41)  are  obtained  by  evaluating  the 
contributions  given  in  Table  (19)  using  the  T intermediates  so 
defined.  The  appropriate  nth-order  contributions  from  the  remaining 
terms  in  Equation  (41)  are  obtained  in  a straightforward  fashion  by 
using  various  low-order  iterates  of  the  unperturbed  t-amplitudes,  the 
X-amplitudes  and  the  derivative  t-amplitudes.  We  summarize  these 
remaining  contributions  for  the  MBPT(2),  MBPT(3),  MBPT(4)  and  D- 
MBPT(“)  models  in  Tables  (20)  through  (23),  respectively.  The  low- 
order  iterates  of  the  derivative  t-amplitudes  are  summarized  in  Table 
(24). 
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Table  19.  The  AO/AO,  AO/MO  and  MO/MO  Derivative  Parts  of  the  Second 
Derivative  of  the  CC/MBPT  Energy. 

aB 

AO/AO  derivative  part  = E <y(a,  (3)v(a,  8)  | p(a,  8)o(a,  8)>  r(yvpa) 

y,  v, p, a 


AO/MO  and  MO/MO  = z r(ij  ab)  <ij  | | ab><txP)  + Z T(ab,  cd)  <ab  | | cd>(<xp) 

derivative  parts  a,b, 

a, b c,d 

+ Z T(ij,kl)  <ij||kl>(aP)  + Z T( ja, bi ) <ja||bi>(ap) 

i > j > j>a» 

k,l  b,i 

+ Z T(ja,bi)  <ja| |bi>(aP)  + E T(ai , be)  <ai | | bc>(a0) 


i»a, 

j.b 

+ E T(jk,ia)  <jk| |ia>(aP) 
i,a, 
j A 


i,a, 

b,c 


where  we  define 

<pq I |rs>(aP)= 

E <y(a,  8)  vi 
y,  v,  p, a 


p(a,  3)a(a,  3)>a 

( (C  C - 
( v pr  as 

C C ) 
ps  ax' 

c _ c U? 

yt  vq  tp 

+ 

(C  C - 
v pr  cs 

C C ) 
ps  ar' 

c c „ uf 

yp  vt  tq 

+ 

(C  _ C 
v pt  as 

C C .) 
ps  at 

C C U? 
yp  vq  tr 

+ 

(C  C - 
v pr  at 

C „ C ) 
pt  ar' 

C C U? 
yp  vq  ts 

+ S ( <tu| |rs>  Ua  UP  + <tq| |us>  UatpUPur  + <tq| |ru>  u“p  U^s 
t,u 

+ <pt||us>  u“q  U*.  + <pt | | ru>  u“q  U^s  + <pq||tu>  u(Js  ) 


+ all  of  the  above  with  a»3 

Note:  The  <yv|  pa>a^  denotes  the  second  derivative  of  an  ordinary  AO 

two-electron  integral,  evaluated  at  a=8=0.  <yv| pa>  is  the 
derivative  with  respect  to  a,  evaluated  at  o=3=0.  r(yvpa)  is 
defined  in  Table  (7). 

See  Table  (6)  for  the  definitions  of  the  T intermediates. 

The  CPHF  coefficients  Ua  are  given  in  Equations  (28,29). 
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Table  20.  Remaining  Contributions  to  the  MBPT(2)  Second  Derivative. 


1 

5 

<ab | | ij> 

“<5<d  * 

l 

5 

<ab||ij>P  tjJ(l) 

1 

+ 2 

*S<i> 

£jk6 

+ 

2 fjkB^J(1) 

1 

' 2 

t“u) 

£cb6 

- 

2 fcba 

Note:  Summation  over  all  labels  is  implied. 
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Table  21.  Remaining  Contributions  to  the  MBPT(3)  Second  Derivative. 


+ 


+ 


+ 


+ 


+ 


\ <ab|  | ij>ex  **b(2)  + \ <ab|  |ij>P  tab(2) 


1 

2 

4tki<2> 

fjk 

S ab 

TT  . i 
ij 

(1) 

+ 

1 

2 

A„ab 

4,ki 

(2) 

fjk 

a ab 

*ij 

(1) 

1 

2 

*S<‘> 

<3 

(2) 

+ 

1 

2 

t^(D 

£jk“ 

A ab 

(2) 

1 

2 

AtJJ(2) 

fcb 

P _ab 

Tij 

(1) 

- 

1 

2 

“3 

(2) 

fcb 

a ab 
. . 
ij 

(1) 

1 

2 

t?*(l) 

ij 

f hP 
cb 

At** 

13 

(2) 

- 

1 

2 

t-O) 

f ' “ 
cb 

A-ab 

(y*.  . 
ij 

(2) 

1 

8 

<cd  | 

|ab>^ 

Tab 

1J 

(1) 

+ 

1 

8 

*3 

(1) 

<cd  | | 

ab>“ 

*ab 

ij 

(1) 

1 

8 

t^d) 

<ij  1 

|kl>P 

(1) 

+ 

1 

8 

«S 

(1) 

<ij  1 1 

kl>“ 

(1) 

ft™ 

<jc| 

|bk>g 

i ab 

Tij 

(1) 

+ 

«s 

(1) 

<jc|  | 

bk>“ 

(1) 

Note:  Summation  over  all  labels  is  implied. 
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Table  22.  Remaining  Contributions  to  the  MBPT(4)  Second  Derivative. 


1 

5 

<ab||ij>a  «®J(3)  + \ <ab||ij>P  tjJ(3) 

+ 

1 

8 

At?j (2)  <cd||ab>P  tjJ(l) 

+ 

1 At^(2)  <cd|  |ab>a  ^(1) 

+ 

1 

8 

t?j(l)  <cd|  |ab>P  At^(2) 

+ 

l tjJ(l)  <cd||ab>a  t*lb.(2) 

+ 

1 

8 

AtjJ(2)  <ij| |kl>P  t^(1) 

+ 

g AtjJ(2)  <ij||kl>a4*J(l) 

+ 

1 

8 

tjJ(l)  <ij||kl>e  At^(2) 

+ 

l <ij||kl>a  A*aJ(2) 

+ 

Atg(2)  <jc|  |bk>P  t»J(l) 

+ 

At£a(2)  <jc| |bk>a  **J<1) 

+ 

tJJ(l)  <jc||bk>P  AxjJ(2) 

+ 

tJJ(l)  <jc|  |bk>a  *ejJ(2) 

+ 

1 

5 

tjJ(l)  Q®j(I,T(l)x^(l)) 

+ 

J tjJ(l)  Q^j (I*^(1)xt(1)) 

+ 

1 

5 

tjJ(l)  QjJ(IP,t(l)xT(l)) 

+ 

\ tjJ(l)  Qj j (IP» T(l)xt(l) ) 

+ 

1 

5 

tjJ(l)  Qjj(Ia,t<l)*e(l)) 

+ 

5 tjJ(l)  Qjj (Ia»J^(l)xt(l)) 

+ 

tj(2)  fbaP  t»(2) 

+ 

t i (2 ) fba“  <(2) 

- 

tj(2)  f..P  tj(2) 

- 

t j (2)  fi:ja  <J(2) 

- 

1 

2 

tjj(l)  <bc | |aj>P  tj(2) 

- 

2 tjj(l)  <bc|  |aj>a  *J(2) 

+ 

1 

2 

tjj(l)  <ib||jk>P  t®(2) 

+ 

2 tjj(l)  <ib||jk>a  **(2) 

+ 

1 

2 

tj(2)  <cj||ab>P  t^(1) 

+ 

2 tj(2)  <cj|  |ab>“  ^(D 

- 

1 

2 

tj(2)  <ij||kb>P  T^d) 

- 

2 tj(2)  <ij||kb>a  **J(1) 

+ 

1 

5 

t£ij<2)  <cd| |bk>3  xjJ(l) 

+ 5 t^®(2)  <cd||bk>a^J(l) 

- 

1 

5 

tkii(2)  <jcMkl>e 

" 5 tJ®J(2)  <jc||kl>a  *»J(1) 

- 

1 

5 

t^(l)  <dk||bc>e  t^(2) 

- 1 tjj(l)  <dk||bc>a^(2) 

+ 

1 

5 

tjJ(l)  <jk||lc>P  Ti j^(2> 

+ l tjJ(l)  <jk!|lc>a^k(2) 

- 

1 

12 

tabcf2.  r 3 xa^c(2) 
tliju;  rkl  Tijku; 

1 abc.0v  f a abc.9. 

- 12  tlij(2>  £kl  *ijk(2) 
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Table  22 — continued 


+ 


1 

12 


dab 

ijk 


<2>  fdc 


3 


abc 

Tijk 


(2) 


dab 

4jk 


(2)  f 


dc 


abc 


(2) 


Note:  Summation  over  all  labels  is  implied. 
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Table  23.  Remaining  Contributions  to  the  LCCD  or  (D-MBPT(®))  Second 
Derivative. 


+ 


+ 


+ 


+ 


\ <ab||ij>a*ab(»)  + 5 <ab||ij>P  tJJ(-) 

1 «.»*>,.  f p rah(^  4-  1 tabf«o  f a *ab(®) 

2 tki(")  fik  Ti1(  > + 2 £ki(  > £jk  *ij^  ; 


«S<-)  £cbS  - 2 *u‘->  fcb“ 


ab. 
jk  ‘ij 

o _ab. 

ij 


1 

2 ‘ijv  ’ ‘cb  ‘ij v 7 2 'ij 

J tjJ(-)  <cd|  |ab>P  xab(»)  + l <cd|  |ab>a 

t^(“)  <ij||kl>P  xab(»)  + 5 t^(“)  <ij||kl>a^b(”) 

tJJ(-)  <jc| |bk>P  xab(»)  + t£a(»)  <jc| |bk>a  *ab(») 


1 

8 lkl 


ij 


Note:  Summation^gver  all  labels  is  implied. 

The  symbol  xa.(°°)  represents  the  converged  LCCD  derivative  t- 
amplitudes.  1'* 
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Table  24.  The  Low-order  Iterates  of  Ta. 


J T^(l)  <ab||cd>  + | X^CD  <ij||kl> 

+ P(ab|ij)  T^(l)  < j c | | bk> 

+ 2 t?j (1)  <ab|  |cd>a  + j tjJ(l)  <ij||kl>a 

+ P(ab|ij)  tg(l)  <jc||bk>a  ] / 


+ P(ab | ij ) AtJJ(2)  <jc| |bk> 

+ j At?j (2)  <ab | | cd>a  + \ At*J(2)  <ij||kl>a 

+ P(ab | ij ) At£J(2)  <jc| |bk>“ 

+ P(ij ) xj(2)  <cj | | ab>  - P(ab)  x®(2)  <ij||kb> 

+ P(ij ) t?(2)  <cj | | ab>a  - P(ab)  tj(2)  <ij||kb>a 
+ \ P(ab)  t^(2)  <cd | | bk>  - \ P(ab)  xgj(2)  <jc||kl> 

♦ \ P(ab)  tj*®<2)  <cd||bk>“  - \ P(ab)  tjjj(2)  <jc||kl>“  ]/ 


M"(3;0)  - t(l)xT(l))  * Qj,(I,  T(l)xt<l)) 

J ^ J 

+ Q*j(Ia,  t(l)xt(l))}/ 


At?j(2)  <ab | | cd>  + \ Ax^(2)  <ij||kl> 
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Table  24 — continued 


[ - \ <bc|  |aj>  + \ xJj(D  <ib|  |jk> 

- \ tb?(l)  <bc | | aj>a  + \ tb*(l)  <ib||jk>a  ] / 


x^(2)  = [-  P(a/bc | k/ij ) xf®(l)  <dk||bc> 


+ P(c/ab | i/jk)  t*b(l)  <jk| |lc> 


- P(a/bc|k/ij)  tfj(l)  <dk||bc>a 

+ P(c/ab | i/jk)  t*b(l)  <jk||lc>a  ] / 


Note:  Summation  over  all  labels  is  implied. 

The  symbol  T represents  the  derivative  t-amplitudes , where  the 
derivative  is  with  respect  to  a,  evaluated  at  a=3=0. 

P(*)  refers  to  various  operators  which  generate  terms  with 
permuted  labels.  See  the  notes  in  Tables  (3)  and  (4)  for  the 
definitions  of  such  operators. 

See  Table  (13)  for  the  low-order  iterates  of  T and  definitions  of 
other  related  quantities. 


CONCLUSION 


The  theory  of  analytic  CC  and  MBPT  first  and  second  derivatives 
has  been  presented;  an  efficient  strategy  has  been  described  for  the 
evaluation  of  molecular  gradients  and  first-order  one-electron 
response  properties.  The  expressions  for  the  first  derivative  of  the 
CCSDT  energy— including  the  CCSDT  response  density— have  been  made 
explicit.  The  expressions  for  the  second  derivative  of  the  CCSDT-1 
energy  have  been  presented  in  terms  of  the  response  density  and 
derivative  t-amplitudes.  The  theory  of  first  and  second  derivatives 
of  the  MBPT  models  through  fourth  order  has  been  extracted  as  a 
special  case  of  the  theory  of  CC  derivatives. 

A report  of  several  analytic  response  properties  of  water  computed 
through  the  SDQ-MBPT(4)  level  of  theory  and  using  a large  Gaussian 
basis  set  has  been  accepted  for  publication.  We  have  published  a 
report  of  the  application  of  MBPT(3)  gradients  to  small  closed-shell 
molecular  systems  using  basis  sets  of  limited  size.  In  the  report, 
we  presented  the  theory  of  MBPT(2)  and  MBPT(3)  response  densities  for 
the  first  time  and  included  values  of  MBPT(2)  and  MBPT(3)  dipole 
moments  as  examples  of  analytic  property  evaluation. 

One-electron  response  properties  through  the  SDQ-MBPT(4)  level  of 
theory  and  for  the  LCCD  model  are  currently  possible.  Molecular 
gradients  are  available  through  the  D-MBPT(4)  model  and  for  the  LCCD 
model.  It  is  possible  to  obtain  plots  of  the  MBPT  response  density, 
as  well;  a recent  study  of  the  dissociation  of  methane  includes  a 
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series  of  MBPT  density  plots  at  various  C-H  bondlengths.  All 
implementation  has  been  general;  that  is,  applicable  to  the  case  of 
UHF  orbitals  as  well  as  RHF  orbitals.  Naturally,  the  UHF-MBPT 
response  density  is  now  being  used  in  spin-density  studies.  Current 
progress  is  being  made  to  extend  gradient/property  evaluation  to  the 
full  MBPT(4)  model  and  the  CCD  model. 
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